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EXECUTIVE  SUMMARV 

This  report  develops  and  describes  the  analytical  models  and  digital 
computer  simulations  of  the  electric  propulsion  systems  for  flywheel 
energy  storage  vehicles.  The  simulation  described  facilitates  a critical 
evaluation  of  candidate  flywheel -powered  electric  vehicles  operating 
over  urban  transit  routes.  The  computer  program  provides  dynamic  or 
time -varying  information  on  voltages,  currents,  tractive  effort,  flywheel 
speed,  etc.  , as  the  vehicle  traverses  a given  roadway  profile  and  con- 
forms to  a specified  velocity  profile.  The  specific  flywheel  electric  - 
drive  propulsion  system  is  well  defined  in  the  computer  model,  and  is 
based  on  manufacturers'  data  for  the  actual  hardware.  The  power  losses 
in  the  major  components  are  computed  at  each  time  interval  and  used  as 
inputs  for  the  thermal  models,  which  determine  the  temperature  rises 
of  the  components.  The  wayside  recharge  of  the  flywheel,  using  the 
on-board  alternator  as  a motor,  is  simulated  and  the  use  of  regenerative 
braking  to  recover  energy  and  store  it  in  the  flywheel  is  also  featured. 

This  simulation  can  be  used  to  determine  whethe r or  not  the  component 
ratings  have  been  properly  selected,  based  on  the  temperature  rises 
observed  after  a run  of  several  round  trips  or  mission  cycles.  Scaling 
laws  are  provided  to  change  the  loss  models  and  the  thermal  models 
in  accordance  with  changes  in  the  machine  ratings  or  sizes. 
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The  simulation  can  also  be  used  to  determine: 

1.  If  the  flywheel  size,  i.  e.  energy  storage  capacity,  is  adequate 
for  the  vehicle/  mission  profile. 

2.  The  energy  efficiency  of  the  propulsion  system  in  terms  of 
average  kWh  required  per  vehicle  route  mile. 

3.  The  effectiveness  of  regenerative  braking  energy  recovery  in 
reducing  the  propulsion  energy  requirement,  i.  e.  is  the  extra 
complexity  of  regenerative  braking  worthwhile  in  terms  of 
energy  savings? 

4.  The  time  required  to  recharge  the  flywheel  at  each  stop,  so 
that  inconveniently  long  recharge  intervals  are  identified. 
Alternate  recharge  strategies  can  then  be  simulated  to  mini- 
mize the  longest  recharge  time.  One  such  scheme  makes  all 
recharge  intervals  of  equal  time. 

Finally,  this  simulation  can  be  used  to  determine  suitability  of  a given 
system  configuration  for  an  intended  mission.  Also,  changes  in  the 
ratings  of  components  and  different  configurations  of  components  can 
be  analyzed  in  order  to  optimize  the  drive  system.  Furthermore,  the 
simulation  method  can  be  applied  to  electric -drive  systems  for  battery 
energy-storage  vehicles,  to  hybrid  vehicles,  and  to  vehicles  operating 


from  wayside  power. 
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1.  INTRODUCTION 


Research  Context 

Development  of  a bus  propelled  by  a flywheel  instead  of  by  an  engine  or 
storage  battery  has  long  interested  urban  transportation  planners. 

The  kinetic  energy  stored  in  a flywheel  supplants  the  need  for  petroleum - 
derived  fuels  and  minimizes  negative  environmental  impacts.  Flywheel- 
driven  vehicles  depend  on  fossil-fuel,  nuclear,  or  hydroelectric  central 
station  power.  They  have  no  emission  fumes  and  low  noise  levels. 

The  flywheel  propulsion  concept  is  not  new.  As  early  as  1935,  an 
Electrogyro  Bus*  powered  by  a kinetic  energy  wheel  was  placed  in 
service  by  Oerlikon  of  Switzerland  (see  Fig.  1).  Designs  that  combine 
flywheel  energy-storage  with  conventional  electric -drive  systems, 
however,  have  not  been  actively  developed  until  recently.  Renewed 
interest  in  this  type  of  vehicle  stems  from  the  dual  concern  with  reduc- 
ing pollution  from  fuel-burning  vehicles  and,  following  the  1974  oil 
embargo,  with  developing  alternative  energy  sources. 

Basically,  a flywheel  is  a spinning  disc  which  stores  kinetic  energy. 
Carried  on  board  a rubber-tired  vehicle,  it  can  be  charged  or  recharged 
using  electric  power  at  selected  wayside  stations.  The  stored 
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Oerlikon  Electrogyro  Bus 


* Parenthetical  names  refer  to  major  components  of  the  baseline  system 
modeled  in  this  report.  While  roughly  equivalent  to  those  used  in  the 
Swiss  bus,  these  components  have  been  varied  somewhat,  for  reasons 
which  will  be  explained  in  Section  2 of  the  report.  (Solid-state  con- 
trollers, for  example,  were  not  available  in  1935. 
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energy  is  expended  between  stations  to  supply:  (1)  traction  power  for 
propelling  the  vehicle  along  the  roadway;  (2)  power  for  overcoming  the 
power  losses  incurred  by  the  major  system  components;  and(3)  auxil- 
iary power  for  on-board  lighting,  heating,  cooling,  and  other  passen- 
ger comforts.  The  sum  of  these  three  power  components  is  the  total 
power  delivered  by  the  flywheel  at  each  instant  of  travel  time.  The 
energy  required  to  propel  the  vehicle  between  charging  stations,  plus 
the  losses  and  auxiliaries,  must  be  equal  to  the  energy  drain  on  the 
flywheel  during  the  same  time  interval.  The  main  purpose  of  the  drive 
system  is  to  convert  the  energy  stored  in  the  flywheel  to  useful  propul- 
sion power  at  the  vehicle  drivewheels. 

Flywheel  energy-storage  systems  are  now  under  serious  consideration 

for  use  not  only  in  buses  but  also  in  subway  cars,  trolleys,  and  other 

vehicles.  For  example,  an  R-32  energy-storage  unit  built  by  the 

Garrett  Corporation  is  currently  being  evaluated  in-service  on  selected 

2 

New  York  City  subways.  Also,  the  Urban  Mass  Transportation  Admin- 
istration recently  sponsored  the  design  of  a flywheel  electric-drive 
system  to  propel  trackless  trolley  coaches  for  the  San  Francisco 

3 

Municipal  Railway.  Hybrid  electric -drive  systems  for  off-highway 

4 

vehicles  have  been  seriously  studies  by  the  U.  S.  Army. 
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Further  development  of  flywheel-driven  vehicles  can  be  facilitated  by 
establishment  of  criteria  and  procedures  for  measuring  the  perfor- 
mance of  alternative  proposed  systems.  Since  the  energy  stored  in 
the  flywheel  must  be  sufficient  for  the  vehicle  to  reach  the  next 
charging  station,  effective  design  procedures  will  be  needed  to  accur- 
ately compute  the  sums  of  the  three  power  and  energy  requirements 
cited  earlier.  In  addition,  procedures  for  selecting  properly  sized 
components  for  the  system  configuration  will  require  realistic  assess- 
ment of  (1)  both  peak  power  delivered  and  maximum  power  dissipated 
as  losses  within  these  components;  and  (2)  the  temperature  rises 
resulting  therefrom. 

Techniques  for  calculating  both  the  required  propulsion  power  and  the 
required  auxiliary  power  of  any  electric -drive  system  are  already  well 
established.  Flywheel  electric-drive  systems  are  more  complex, 
however,  than  the  conventional  electric-drive  systems  used  for  battery- 
and  wayside -powered  vehicles.  Application  of  standard  design  techniques 
could  result  in  erroneous  estimation  of  (1)  flywheel  energy-storage 
requirements  and  (2)  the  temperature  rises  of  the  system  components. 
Despite  the  cited  advantages  of  flywheel  propulsion,  the  losses  in  the 
major  components  of  such  systems  often  constitute  a substantial  portion 
of  the  output.  Procedures  more  accurate  than  the  standard  design 
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techniques  will  thus  be  needed,  first,  to  calculate  these  losses  and 
then,  to  use  these  calculations  for  determining  the  total  energy 
required  from  the  flywheel  and  the  temperature  rises  of  the  major 
system  components. 

Performance  Criteria  for  Flywheel  Electric-Drive  Systems 

This  report  can  be  used  for  evaluating  the  degree  to  which  the  config- 
uration of  major  components  and  the  component  ratings  selected  for 

5 

a proposed  propulsion  system  will  optimize  its  performance.  The 
report  presents  a procedure  for  developing  computer  models  of  the 
electric -drive  system  for  a flywheel  energy-storage  vehicle.  The 
computer  simulates  the  overall  operation  of  the  vehicle  and  its  drive 
system  as  well  as  the  electromechanical  and  thermal  operation  of  the 
drive  system  components.  The  "baseline"  system  modeled  was  pat- 
terned after  an  existing  vehicle  traveling  a known  route  in  accordance 
with  a predetermined  velocity  profile.  Any  parameter  of  the  route,  of 
the  vehicle,  or  of  the  drive  system  components  could,  however,  be  varied 
and  the  overall  effect  on  system  performance  easily  determined.  These 
modeling  techniques  can  be  applied  to  conventional  battery-powered 
as  well  as  flywheel  electric-drive  systems,  and  to  any  vehicle  operating 


from  wayside  power. 
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To  evaluate  the  performance  of  a proposed  system,  two  primary  factors 
should  be  considered: 

Proper  ratings  of  the  major  components 

Adequacy  of  the  flywheel  energy-storage  capacity  . 

The  ratings  selected  for  the  major  components  determine  how  heavily 
they  are  loaded  in  actual  use.  If  a selected  rating  is  too  small,  the 
undersized  unit  is  overloaded  and  operates  less  efficiently,  i.  e.  , with 
larger  power  losses.  This  larger  loss  is  dissipated  in  the  form  of 
heat  within  the  component  and  thus  produces  a higher  temperature  rise 
than  "allowable  by  the  manufacturer.  " The  excess  temperature  shortens 
the  life  expectancy  of  the  unit.  * Allowable  temperature  rises,  deter- 
mined on  the  basis  of  the  class  of  insulation  used,  typically  range  from 
80  C to  160°  C.  As  a rough  rule-of-thumb,  for  every  10°  C a component 
operates  over  the  allowable  temperature,  its  life  expectancy  is  halved. 

If  a selected  rating  is  too  large,  the  temperature  rise  is  minimal  but 
the  penalties  of  increased  size,  weight,  and  cost  reduce  the  attractiveness 
of  the  design.  Optimal  sizing  of  a major  component  is  achieved  when  the 
final  temperature  rise  that  occurs  after  several  round  trips  is  5 to  10 


* Sudden  and  complete  failures  of  major  components  are  also  possible 
when  the  allowable  temperatures  are  exceeded.  The  likelihood  of 
such  breakdowns  depends  on  the  type  of  major  component  in  question 
(e.  g.  motor  or  solid-state  device)  and  the  degree  of  excess  tempera- 
ture. 
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degrees  below  the  manufacturer's  allowable  temperature  rise. 
Misspecification  of  the  flywheel  size  can  have  varied  effects:  under- 
specifying  will  not  allow  for  sufficient  energy  to  propel  the  vehicle 
from  station  to  station;  overspecifying  will  increase  flywheel  weight 
at  the  expense  of  the  vehicle's  passenger  capacity.  Accurate  deter- 
mination of  the  flywheel  size  in  the  early  stages  of  design  is  particu- 
larly important  since  correcting  a misspecification  involves  time- 
consuming  and  costly  changes  in  vehicle  design. 

This  time  and  cost  constraint  is  considerably  more  stringent  than  that 
on  the  design  of  other  forms  of  propulsion.  For  example,  minor  mod- 
ifications can  usually  accomplish  major  changes  in  the  power  ratings 
of  batteries  or  engines.  If  underpowered,  a vehicle  that  operates  on 
electric  batteries  can  have  more  batteries  added;  one  powered  by  a 
gasoline  engine  can  have  minor  changes  made  to  the  spark  advance  or 
air  intake.  Any  flywheel  modifications,  by  contrast,  would  entail  not 
only  changing  the  containment  ring  used  to  isolate  the  flywheel  in  case 
of  mishap  (and  thus  the  size  of  the  envelope)  but  also  redesigning  the 
vehicle  layout. 

Optimal  flywheel  sizing  procedures  should  determine  the  energy- 
storage  capacity  needed  by  accurately  computing  all  known  energy 
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requirements.  A margin  of  40  to  50  percent  should  then  be  added 
to  these  requirements  to  allow  for  contingencies,  such  as  high  winds, 
and  to  prevent  flywheel  spin-down  below  half  speed. 

The  key  to  designing  a system  that  fulfills  all  performance  criteria 
is  careful  determination  of  the  power  losses  of  each  major  system 
component.  Power  losses  directly  determine  the  temperature  rises 
of  the  components;  the  losses  are  required  to  derive  the  component 
efficiencies  and  to  calculate  the  energy  drain  from  the  flywheel. 
Inaccurate  determinations  of  losses,  then,  can  lead  to  poor  sizing 
both  of  the  major  components  and  of  the  flywheel.  In  the  computer 
simulation,  power  losses  are  calculated  at  each  instant  of  time,  as  a 
function  of  the  dynamic  state  of  the  system. 

Two  additional  criteria  evaluated  by  the  simulation  are:  usage  con- 
straints, particularly  the  time  required  to  recharge  the  flywheel  at  a 
station;  and  the  overall  efficiency  by  which  flywheel  energy  is  converted 
to  useful  propulsion  power.  Energy-conversion  efficiency  is  clearly 
important  in  a time  of  increasing  costs  for  all  forms  of  energy  and 
will  influence  the  projected  operating  costs  for  the  vehicle.  Although 
both  these  factors  are  usually  determined  only  from  actual  test  data. 
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computer  evaluation  is  possible  primarily  as  a byproduct  of  the 
"dynamic"  simulation  described  in  the  next  section. 

Simulation  Approach 

The  main  computer  program  is  divided  into  two  major  portions. 

The  first  portion  simulates  the  dynamic  behavior  of  the  vehicle  as  it 
makes  a run  over  the  route.  Using  as  input  the  velocity  profile  of 
the  vehicle  and  the  grade  profile  of  the  roadway,  the  computer  calcu- 
lates at  each  time  increment:  the  electrical  and  mechanical  quantities 
of  the  drive  system  (such  as  power,  voltage,  current,  and  speed);  the 
required  propulsion  power;  the  power  losses  of  the  major  components; 
and  the  power  required  from  the  flywheel.  The  calculated  power  losses 
are  used  to  determine  the  efficiencies  of  the  components  at  any  given 
instant;  flywheel  power  requirements  are  summed  while  the  vehicle  is 
traveling  to  find  the  flywheel  energy  drain  (spin-down)  between  stations. 

In  the  second  thermal  portion  of  the  main  program,  the  power  losses 
already  determined  are  used  as  inputs  to  thermal  models  in  order  to 
compute  the  temperature  rises  of  individual  components.  Thermal 
models  are  developed  only  for  those  major  components  whose  tempera- 
ture rises  are  expected  to  be  significant. 
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These  modeling  procedures  are  more  applicable  to  the  design  of 
flywheel  energy-storage  vehicles  than  techniques  traditionally  used 
for  electric-drive  systems.  The  critical  data  requirement  for  design- 
ing a properly-sized  flywheel  is  the  total  energy  required  from  the 
flywheel.  The  methods  used  in  the  Flywheel  Propulsion  Simulation  to 
calculate  this  quantity  as  well  as  to  evaluate  the  selected  ratings  of 
other  major  components  either  expand  on  or  depart  from  the  tradi- 
tional approach  in  a number  of  important  respects. 

Any  approach  to  designing  an  electric-drive  system  starts  with  a 
vehicle  propelled  over  a roadway  in  accordance  with  the  desired  velo- 
city profile.  Data  on  either  an  existing  or  a hypothetical  vehicle  may 
be  used,  and  aerodynamic  drag  may  be  taken  into  account  if  desired. 

In  a first  step,  the  required  tractive  force  is  determined  on  the  basis 
of  these  data  as  a function  of  time,  usually  by  using  a computer  simu- 
lation to  insure  greater  accuracy.  The  required  propulsion  power  is 
also  calculated  as  the  product  of  tractive  force  and  velocity. 

Following  this  step,  however,  the  basic  assumptions  and  calculation 
procedures  of  the  Flywheel  Propulsion  Simulation  differ  significantly 
from  the  traditional  approach.  Traditionally,  component  ratings  are 
selected  by  approximating  both  the  average  and  the  peak  power  that 
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the  component  will  have  to  deliver  in  order  to  get  the  vehicle  over 
the  route  on  schedule.  To  determine  the  average  loading  of  the  electric  - 
drive-system  components,  either  the  time-averaged  or  the  root-mean- 
square  (rms)  value  of  the  propulsion  power  is  used  along  with  an  effi- 
ciency factor,  which  is  assigned  to  each  component  on  the  basis  of 
manufacturers'  test  data  obtained  under  steady-state  rated  conditions. 

This  procedure  assumes  that  efficiency  values  will  remain  meaningful 
as  the  load  undergoes  large  variations  during  the  mission.  To  deter- 
mine whether  the  thermal  heating  of  a component  will  be  within  allow- 
able limits,  it  furthermore  relies  heavily  on  the  designer's  judgment 
of  the  effect  of  average  and  peak  power  losses  on  service  life  of  the 
component.  The  traditional  methodology,  in  sum,  does  not  take  into 
account  the  effects  of  wide  variation  of  operating  conditions  and  losses 
during  a typical  run  of  a vehicle. 

In  contrast,  one  of  the  basic  assumptions  of  this  report  is  that  com- 
ponent efficiency  will  vary  widely  under  the  dynamic  conditions  of  the 
roadway.  Because  the  modeling  techniques  compute  power  losses  at 
each  instant  of  time,  the  Flywheel  Propulsion  Simulation  can  determine 
average  and  peak  losses  more  realistically  and  can  more  accurately 
calculate  on  the  basis  of  these  losses  both  the  flywheel  energy-storage 
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capacity  required  and  the  thermal  heating  that  will  occur  in  the 
major  components. 

The  inclusion  of  detailed  loss  models  for  each  component  and  the 
development  of  thermal  models  for  the  major  components,  thus, 
constitute  a significant  expansion  of  the  technology  available  for 
designing  electric -drive  systems.  Because  the  simulation  includes 
thermal  models,  it  can  both  track  temperatures  and  evaluate  the 
suitability  of  components  of  different  ratings  and  designs  for  inclusion 
in  a proposed  system.  Because  it  includes  highly  detailed  loss  models, 
the  simulation  is  particularly  applicable  to  flywheel  energy-storage 
systems,  where  losses  are  so  significant  a part  of  the  output.  Though 
traditional  techniques  necessarily  take  losses  into  account,  the  actual 
loss  calculations  for  each  component  are  embodied  in  only  one,  approx- 
imated, efficiency  factor. 

It  should  be  noted  that  while  the  accuracy  of  the  modeling  procedure 
derives  from  the  dynamic  approach,  the  credibility  of  the  loss  models 
depends  not  only  on  their  detail  but  also  on  their  use  of  manufacturers' 
loss  data  for  existing  hardware.  The  manufacturers'  efficiency  data;, 
traditionally  employed,  were  not  used  in  this  simulation  except  to 
confirm  the  accuracy  of  the  loss  components  of  the  models. 
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Scope  of  Analysis 

The  Flywheel  Propulsion  Simulation,  in  summary,  can  be  considered 
both  an  evaluation  tool  and  a design  tool.  On  one  hand,  it  can  be  used 
for  comparing  proposed  flywheel  electric-drive  systems  with  respect 
to  specific  performance  criteria.  On  the  other  hand,  it  can  deter- 
mine selection  of  the  correct  ratings  for  major  system  components. 

Such  changes  in  the  size  of  units  do  not  necessarily  require  basic 
redesign  of  the  system,  however.  More  fundamental  changes  in  a 
system  configuration  --  e.  g.  , substituting  an  ac  for  a dc  traction 
motor  or  using  a chopper  instead  of  a phase -controlled  rectifier  -- 
are  also  possible  with  this  simulation  method.  Variations  in  sizes  of 
components  are  therefore  considered,  but  detailed  evaluation  of  alter- 
natives at  this  time,  however,  would  have  demanded  extensive  repro- 
gramming of  the  simulation,  a change  clearly  not  permitted  by  the 
time  constraints  of  this  study. 

Also  excluded  from  the  scope  of  analysis  was  the  possibility  of 
coupling  a flywheel  energy-storage  unit  with  a mechanical-  rather  than 
an  electric -drive  system.  The  baseline  electric-drive  system  modeled 
offered  a number  of  advantages:  the  flexibility  and  high  efficiency  of 
solid-state  control;  the  availability  of  electric  machines  of  reasonably 
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high  efficiencies;  and  greater  flexibility  in  the  arrangement  of  the 
major  components  within  the  vehicle.  It  furthermore  bypassed 
several  disadvantages.  A mechanical-drive  system  would  have 
required  an  additional  auxiliary  power  generator  and  entailed  a more 
complex  and  expensive  process  for  recharging  the  flywheel. 

A description  of  the  baseline  system  modeled  --  the  components  and 
prototype  existing  vehicle  as  well  as  the  modes  of  system  operation  -- 
is  presented  in  the  next  chapter  of  this  report,  followed  by  a descrip- 
tion of  modeling  procedures  in  Section  3.  The  parameters  for  the 
power  loss  and  thermal  models  are  detailed  in  Section  4.  Because 
the  simulation  results  were  necessarily  constrained  by  the  baseline 
parameters  modeled,  selected  variations  in  those  parameters  were 
considered,  though  not  actually  modeled,  and  are  also  discussed  in 
Section  4.  Results  of  a sample  run  for  the  baseline  system,  as 
well  as  recommendations  for  the  use  of  the  simulation,  and  the 
design  of  electric -drive  systems,  are  given  in  Section  5.  Several 
appendices  provide  the  mathematical  equations  which  are  the  dynamic, 
loss  and  thermal  models  of  the  drive  system;  definitions  of  the  vari- 
ables used  in  the  report;  typical  printouts  of  a computer  run,  and 
manufacturers'  data  on  the  components  for  the  baseline  system. 
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2.  BASELINE  SYSTEM:  COMPONENTS  AND  MODES  OF  OPERATION 


Configuration  of  Major  Components 

The  baseline  flywheel  energy-storage  vehicle  was  patterned  after  an 
existing  vehicle:  the  21 -passenger,  rubber-tired  PRT  (personal  rapid 
transit  system) ^ vehicle  built  by  Boeing  for  the  Morgantown,  West 
Virginia  transportation  project,  which  was  sponsored  by  the  Urban  Mass 
Transportation  Administration.  This  vehicle  (see  Fig.  2)  was  selected  for 
modeling  mainly  because  it  is  already  operant.  While  the  propulsion 
system  modeled  for  the  baseline  vehicle  has  several  components  in 
common  with  the  existing  vehicle,  its  configuration  is  more  complex 
since  it  adds  a flywheel  and  an  alternator. 

The  configuration  of  components  selected  to  develop  the  modeling 

procedure  is  that  most  frequently  proposed  by  other  designers  for 

7 

flywheel  energy-storage  vehicles.  It  includes  five  major  units  (see 
Fig.  3),  which  serve  the  following  general  functions: 

The  flywheel  stores  energy  that  is  charged  into  it  at  a 
wayside  power  connection  and  delivers  energy  when  the 
vehicle  is  traveling  between  stations.  The  component 
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Fig.  2 

Morgantown  Vehicle 
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modeled  is  characterized  by  only  two  parameters: 

maximum  speed  (12,000  r/  min),  and  moment  of  inertia 
2 

(45  slug-ft  ),  equivalent  to  13.4  kWh. 

The  alternator  directly  coupled  to  the  flywheel, 
supplies  the  auxiliary  load  and  generates  ac  power  for 
vehicle  propulsion  as  the  flywheel  slows  down.  During 
recharge,  it  accepts  ac  wayside  power  as  a motor. 

The  component  modeled  is  a Westinghouse  brushless, 
3-phase,  4-pole,  synchronous  machine  that  is  rated  at 

g 

75.8  kVA  unity  power  factor  and  12,000  r/min.  The 

alternator  was  selected  to  match  the  power  rating  of  the 

dc  traction  motor  and  the  maximum  speed  of  the  flywheel. 

The  alternator-flywheel  is  similar  to  that  for  the  San 

Francisco  Municipal  Railway  prototype  flywheel  energy- 

3 

storage  trolley  coach. 

The  solid-state  controller  consists  of  a 6-thyristor  (SCR), 

3 -phase,  full-wave  phase-controlled  bridge  rectifier. 

It  normally  modulates  and  rectifies  the  ac  power  for  use 

in  the  dc  motor  but  may  also,  during  braking,  convert 

9 

dc  to  ac  power.  The  F-400  Power  Pack  manufactured 
by  Power  Semiconductors,  Inc.  was  modeled. 
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Major  Components  of  Baseline  Propulsion  System 
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The  dc  traction  motor  provides  the  torque  to  drive  the 
vehicle.  The  component  modeled  is  a compound-wound, 
4-pole  machine  that  is  rated  at  70  hp  and  2,  730  r/  min. 
This  type  ASEA-built  motor  was  installed  by  Rand- 
tronics*^  in  the  Morgantown  vehicle. 

The  single-stage  transmission,  a simple  gearbox, 
provides  a speed  reduction  of  7.  17:1  as  it  couples  the 
dc  motor  to  the  vehicle  drive -wheels. 

The  first  two  components  comprise  the  energy-storage  unit;  the  last 
three  can  be  called  the  "power  train,  " which  is  subdivided  into  elec- 
trical and  mechanical  parts.  The  transmission  and  drive-wheels 
comprise  the  mechanical  part.  As  the  vehicle  is  propelled  along  the 
roadway,  the  direction  of  power  is  usually  from  the  energy-storage 
unit  via  the  power  train  to  the  drive -wheels.  During  braking,  if  that 
operation  is  performed  electrically  rather  than  with  friction  brakes, 
this  direction  of  power  is  reversed. 

The  data  used  in  the  computer  simulation  for  these  components  were 
derived  from  two  sources:  for  the  dc  traction  motor  and  gearbox, 
from  published  data  on  the  components  actually  used  in  the  Morgan- 
town vehicle;  for  the  alternator  and  solid-state  controller,  from 
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manufacturers'  equipment  catalogs  and  from  personal  communica- 
tions with  manufacturers'  personnel.  Manufacturers'  data  on  all 
components  are  detailed  in  Appendix  E. 

Of  the  components  modeled,  the  ac  and  the  dc  machines  typically 
have  the  lowest  efficiencies  --  each  about  0.  88  at  base  conditions 
and  considerably  lower  under  adverse  combinations  of  load  and 
speed.  Efficiency  of  the  gearbox  can  be  expected  to  remain  above 
0.  90  at  all  times;  of  the  solid-state  controller,  above  0.  97.  Hence, 
the  two  components  with  lower  efficiencies  (i.  e.  , larger  power 
losses)  required  more  elaborate  loss  models;  the  loss  models  for 
the  gearbox  and  the  solid-state  controller  were  much  less  complex. 
The  flywheel  losses  were  neglected;  they  could  be  included  as  part 
of  the  alternator  friction  and  windage  losses. 

The  Mission  Profile 

The  Morgantown  PRT  vehicle  modeled  was  assumed  to  traverse  the 
route  between  the  Walnut  and  Engineering  stations  on  which  that 
vehicle  is  presently  operating.  The  Flywheel  Propulsion  Simulation 
adopted  both  the  route  elevation  and  the  vehicle  velocity  profiles 
currently  in  use  (see  Fig.  4).  The  vehicle  covers  a distance  of  2.  065 
mi  one-way  (4.  13  mi  round-trip)  and  travels  at  a maximum  speed  of 
30  mi/  h.  One-way  travel  time  between  the  two  stations,  which  are 


Elevation  Profile 

Walnut  Station  Engineering  Station 

(Point  A)  (Point  B)  , 
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Roadway  Elevation  and  Vehicle  Velocity  Profiles  for  the  Morgantown  Route 
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labeled  A and  B respectively,  is  roughly  seven  minutes.  The  portion 
A-to-B  is  the  outbound  segment  of  the  route;  B-to-A,  the  inbound  segment. 

The  complete  round-trip  between  these  stations  can  be  called  the 
"mission  cycle"  of  the  vehicle.  Initially,  the  flywheel  is  charged  to 
maximum  speed  at  the  Walnut  station  (Point  A).  The  vehicle  then 
travels  under  flywheel  propulsion  to  the  Engineering  station  (Point  B), 
where  the  flywheel  is  recharged  during  passenger  boarding.  Finally, 
the  vehicle  returns  to  the  starting  point,  and  the  mission  cycle  is 
repeated  following  a flywheel  recharge  at  Point  A. 

In  the  computer  simulation,  both  mission  profile  parameters  and 
vehicle  parameters  are  used  to  compute  the  required  tractive  force 
and  the  required  propulsion  power  at  each  instant  of  travel.  The  three 
mission  profile  parameters  are: 

Roadway  grade  (slope  of  roadway  elevation  plot) 

Vehicle  velocity 
Wind  velocity  . 

Prestored  values  of  these  parameters  are  used  at  each  instant  of 
elapsed  travel  time.  The  following  vehicle  parameters  are  prestored 
at  the  start  of  a computer  run: 


Effective  weight 
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Tire  friction  coefficients  (rolling  and  coulomb  friction) 
Aerodynamic  drag  forces  (drag  coefficient  and  frontal  area). 


Headwinds  and/  or  tailwinds  may  be  used  in  the  drag  calculation. 

The  Flywheel  Propulsion  Simulation  is  capable  of  modeling  a more 
elaborate  mission  cycle  (one  with  three  or  more  stations,  for  example) 
and  any  number  of  successive  round-trips  (i.  e.  , extended  missions). 

In  modeling  the  continuous  performance  of  the  vehicle,  this  multiple - 
trip  capability  is  needed  to  find  the  final  temperature  rises  of  those 
components,  i.  e.  , the  electric  machines  that  have  long  thermal  time 
constants. 

Modes  of  Operation 

The  flywheel  propulsion  system  operates  alternately  in  two  different 
modes.  In  the  traveling  mode,  flywheel  energy  is  expended  between 
stations  to  supply  the  three  major  requirements  of  the  system.  In 
the  wayside  charging  mode,  that  energy  is  replaced  at  each  stop  as 
the  flywheel  is  recharged  to  maximum  speed. 

Wayside  Charging  Mode.  The  charging  or  recharging  process  (see 
Fig.  5)  begins  when  the  vehicle  is  parked  at  a station  and  passengers 


Wayside 
Power  Source 
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Fig,  5 

Wayside  Charging  Mode;  System  Diagram 
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are  boarding.  The  dc  traction  motor  and  solid-state  controller  are 
electrically  disconnected;  the  auxiliary  load  is  supplied  by  wayside 
power.  A stationary  charging-station  inverter  provides  3 -phase 
adjustable -frequency  electric  power  to  the  alternator  operating  as 
a motor  to  accelerate  the  flywheel. 

As  the  flywheel  speeds  up,  the  input  power  delivered  to  the  ac  machine 
is  controlled  to  be  proportional  to  the  alternator-flywheel  speed  (i.  e.  , 
the  ac  frequency),  while  the  armature  current  and  the  power  factor 
are  held  constant  (see  Fig.  6).  Input  voltage  increases  in  proportion 
to  the  speed  of  the  alternator;  at  maximum  speed,  this  voltage  is 
assumed  to  be  20  percent  higher  than  the  rated  terminal  voltage  of  the 
alternator.  The  20  percent  higher  charging  voltage  is  used  to  reduce 
the  armature  current  and  to  minimize  alternator  overheating  during 
each  recharge. 

The  charging  process  usually  continues  until  the  maximum  (rated) 
speed  of  the  alternator  and  flywheel  ( 1 2 , 000  r/  min)  has  been  reached. 

It  could,  however,  be  stopped  at  any  point  if  additional  constraints  on 
charging  are  introduced  (such  as  an  arbitrary  upper  bound  on  the 
time  available  for  recharging  the  flywheel  at  each  station).  The 
charging  terminated,  the  vehicle  proceeds  down  the  roadway  in  the 


traveling  mode. 


Input  Voltage 
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Alternator  Speed 


* Alternator  operating  as  a motor  at  constant  armature  current 
of  2.  2 pu  and  constant  power  factor  of  -1. 

Fig.  6 


Wayside  Charging  Mode  for  Baseline  System: 


Input  Power  and  Input  Voltage  vs.  Alternator  Speed* 
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A critical  problem  with  flywheel  electric -propuls ion  systems  is  the 
long  recharge  time.  For  the  baseline  system  modeled,  recharge  to 
full  flywheel  speed  and  energy  required  215  s after  the  vehicle  trav- 
eled an  uphill  outbound  segment  and  80  s after  a downhill  inbound 
segment.  The  longer  recharge  interval  is  about  50  percent  of  the 
one-way  travel  time;  the  shorter  interval,  about  20  percent.  How- 
ever, the  charging  intervals  could  be  made  equal,  at  about  157  s 
each,  corresponding  to  39  percent  of  the  one-way  travel  time.  The 
flywheel  would  not  be  fully  charged  when  the  vehicle  starts  its  down- 
hill inbound  segment.  Since  the  alternator  would  be  running  at  less 
than  its  full  speed,  as  well,  its  extra  field-circuit  losses  would  cause 
the  overall  energy  losses  for  that  segment  to  be  greater  than  for  the 
baseline  system.  These  recharge  times  would  probably  be  an  unac- 
ceptable inconvenience  to  the  general  public.  In  fact,  long  recharge 
time  caused  abandonment  of  the  early  Swiss  Electrogyro  Bus.  The 
charging  time  can  be  reduced  by  increasing  the  size  of  the  alternator 
and  by  raising  the  power  delivered  from  the  wayside  to  the  alternator 
and  flywheel. 

Traveling  Mode.  Using  flywheel  energy,  the  vehicle  moves  along  the 
roadway  in  accordance  with  the  predetermined  velocity  profile.  To 
supply  the  correct  tractive  force  required  for  maintaining  this  velocity, 
the  firing  angle  of  the  solid-state  controller  is  adjusted  continuously. 


which  in  turn  determines  the  load  on  the  alternator  at  each  instant  of 
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travel.  These  adjustments  are  required  to  account  for  time -varying 
changes  in  required  velocity  and  tractive  force  that  occur  as  the 
vehicle  travels  between  stations.  For  a given  vehicle  configuration, 
these  variations  result  from  changes  in: 

Inertia  force  (proportional  to  vehicle  acceleration) 

Gravity  force  (proportional  to  roadway  grade  or  slope) 
Tire  friction  force  (nearly  constant-changes  slightly  with 
vehicle  velocity) 

Aerodynamic  drag  force  (proportional  to  square  of  sum 
of  vehicle  velocity  and  wind  velocity). 

In  the  computer  simulation,  then,  the  velocity,  grade,  and  wind 
profiles  are  pre-stored  as  the  mission  profile.  Tractive  force  is 
calculated  at  each  instant  of  travel  as  the  sum  of  four  components: 

1)  inertia;  2)  gravity;  3)  tire  friction;  and  4)  aerodynamic  drag. 

The  field  current  of  the  alternator  is  regulated  to  maintain  constant 
terminal  voltage,  even  though  the  alternator  slows  down  from  12,000 
r/min  to  as  low  as  6,000  r/ min.  The  auxiliary  load  is  thus  supplied 
at  constant  ac  voltage,  but  at  a frequency  from  400  Hz  to  as  low  as 
200  Hz. 

The  traveling  mode  has  two  mutually  exclusive  forms,  either  motoring 
or  braking,  depending  on  the  positive  or  negative  sign  of  the  tractive 
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force  required  at  the  drive -wheels.  During  motoring,  positive  tractive 
force  is  required  to  push  the  vehicle  in  the  direction  of  travel.  When 
the  vehicle  slows  down,  e.  g.  , before  a station  or  on  a steep  downgrade, 
negative  tractive  force  is  required. 

Regenerative  braking,  rather  than  either  dynamic  braking  or 
friction  braking  used  on  the  Morgantown  vehicle  was  selected  for  the 
simulation,  since  it  is  potentially  more  energy-efficient.  The  type 
of  braking  used  in  a flywheel  energy-storage  vehicle  can  influence 
total  energy  consumption  and  therefore  the  amount  of  flywheel  energy 
expended.  Regenerative  braking  can  transfer  some  of  the  vehicle's 
kinetic  energy  to  the  flywheel;  in  both  friction  and  dynamic  braking, 
that  energy  would  be  dissipated  in  the  form  of  heat.  The  energy 
savings  from  regenerative  braking  are  potentially  greatest  for  heavy 
vehicles  (which  have  greater  kinetic  energy)  and  for  vehicles  that: 
stop  frequently  for  passengers  (i.  e.  , over  five  times  per  mile);  require 
high  braking  decelerations;  or  travel  routes  with  steep  grades.  These 
are  conditions  for  the  largest  negative  tractive  force  and  generally  the 
largest  negative  propulsion  power  at  the  drive-wheels- 

During  regenerative  braking,  the  direction  of  power  is  the  opposite  of 
that  during  motoring  (see  Fig.  7).  The  direction  of  the  current  in  the 
field  windings  is  also  reversed,  but  that  in  the  armature  winding 
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Power 

Alternator  Train  Drive  Wheel 


Power 

Alternator  Train  Drive  Wheel 


Regenerative  Braking 
Fig.  7 


Traveling  Mode:  Direction  of  Power  Flow  During 


Motoring  and  Full  Regenerative  Braking 
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remains  the  same  as  in  motoring.  In  motoring,  the  dc  traction  motor 
receives  energy  for  propelling  the  vehicle  from  the  solid-state 
controller;  in  braking,  it  operates  rather  as  a dc  generator,  receiving 
energy  from  the  drive-wheels  and  sending  it  into  the  controller,  where 

it  is  converted  into  ac  power.  This  ac  power  is  used  primarily  to 
provide  the  auxiliary  load.  Any  surplus  power,  after  the  component 

losses  have  been  accounted  for,  is  used  to  drive  the  alternator  as  a 

motor  and  thus,  to  supply  the  torque  that  accelerates  the  flywheel 
during  regenerative  braking. 

Two  types  of  regenerative  braking  --  either  full  or  partial  --  are 
possible.  If  the  traction  power  is  greater  than  or  equal  to  the  auxiliary 
load  plus  the  sum  of  all  the  losses  in  the  major  components,  full  regen- 
erative braking  occurs;  if  not,  the  braking  is  only  partial.  During 
partial  braking,  the  dc  motor  operates  as  a generator;  the  alternator, 
however,  does  not  operate  as  a motor.  Rather,  it  also  operates  as  a 
generator  but  supplies  a reduced  electrical  load.  This  reduction  of 
the  electrical  load  on  the  alternator  diminishes  the  energy  drain  on 
the  flywheel.  Even  in  partial  regenerative  braking,  this  energy  drain 
is  less  than  in  friction  braking. 
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In  sum,  the  direction  of  power  can  be  either  from  the  flywheel  to  the 
drive -wheels  (during  motoring),  or  in  the  reverse  direction  (during 
regenerative  braking).  Regardless  of  the  direction,  the  losses  of 
the  major  components  must  always  be  computed  and  their  effects  on 
flywheel  speed  properly  accounted  for.  In  the  Flywheel  Propulsion 

Simulation,  however,  the  sequence  in  which  losses  are  computed  always 
proceeds  from  the  drive-wheels  to  the  flywheel.  To  determine  the 
power  required  from  the  flywheel,  which  is  also  the  rate  of  flywheel 
energy  discharge,  the  required  traction  power  is  worked  back 
through  the  loss  models.  This  sequence  of  computations  holds  for 
both  the  motoring  and  regenerative-braking  modes;  in  braking, 
however,  the  tractive  force  and  the  traction  power  are  negative. 
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3.  COMPUTER  SIMULATION  PROCEDURES 


Scope  of  Simulation 

The  Flywheel  Propulsion  Simulation  has  been  implemented  in  the 
Digital  Equipment  Corporation  DEC  System -10  computer  located 
at  the  Transportation  System  Center,  Cambridge,  Massachusetts. 
Because  the  simulation  is  programmed  in  Fortran  IV,  a universal 
language,  the  program  can  be  used  on  most  other  large-scale 
machines  with  little  or  no  modification.  This  section  describes  the 
operation  of  the  computer  program  and  the  procedures  for  calcula- 
ting propulsion  power  and  auxiliary  power  requirements.  More 
complex  portions  of  the  simulation,  including  the  power  loss  and 
thermal  models  as  well  as  alternative  ratings  of  major  components, 
are  detailed  in  the  next  chapter;  a sample  run  of  the  program  for  the 
baseline  system  is  described  in  Section  5. 

Operation  of  Computer  Program 

The  Flywheel  Propulsion  Simulation  consists  of  an  executive  program, 
a main  program,  and  nine  subroutines.  The  executive  program  is 
called  VSPC  (Vehicle  System  Performance  Calculator);  it  serves  as 
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the  program  controller,  since  it  calls  in  the  main  program  and 
subroutines  as  needed.  The  main  program  is  called  FLY  (Flywheel); 
it  performs  all  the  computations  required  for  simulating  operation  of 
the  baseline  system,  including  the  following: 

Required  tractive  force  and  propulsion  power 

Power  losses  in  all  major  components 

Total  power  at  flywheel  during  discharge  and  recharge 

Flywheel  speed  and  its  kinetic  energy 

Electric  currents  and  voltages 

Thermal  heating  of  select  major  components  . 

Nine  auxiliary  subroutines  are  used  in  conjunction  with  the  main 
program:  the  two  subroutines  that  process  input  data  are: 

FLY$  1 a listing  of  all  input  constants 
PROF  1 a subroutine  to  process  "raw"  mission  profile 
data  . 

The  five  subroutines  that  collect,  label  or  summarize  the  output  data 
are: 

FORM  1 prepares  comment  statements  for  the  tabular 
dynamic  and  power  loss  printout  (FOR  <4  3) 


FORM  2 


prepares  page  titles  for  the  FOR<i  3 printout 
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FORM  3 prepares  a summary  output  to  follow 
tabular  outputs 

FORM  4 prepares  page  titles  for  the  tabular  thermal 
output  (FOR  26.  DAT) 

FORM  5 prepares  comment  statements  for  the 
FOR  26.  DAT  output  . 

The  tabular  outputs  for  a run  are  stored  in  the  two  subroutines  used 
as  data  storage  files: 

FOR  (j)  3.  DAT  dynamics  and  power  loss  tabular  output 
FOR  26.  DAT  thermal  tabular  output  . 

The  flowchart  for  the  Flywheel  Propulsion  Simulation,  shown  in  Fig.  8, 
indicates  the  overall  sequence  of  computations  and  the  decisions  made 
at  each  time  increment  of  the  simulation.  A computation  step-size  of 
one  second  was  used. 

The  flowchart  contains  two  major  loops,  for  the  wayside  charging  and 

\ 

traveling  modes.  The  program  is  begun  in  the  wayside  charging  loop, 
with  the  flywheel  initially  at  low  speed,  and  remains  there  until  the 
flywheel  reaches  its  rated  or  full  speed.  It  then  switches  to  the  travel- 
ing loop  and  calls  in  new  mission  profile  inputs  (i.  e.  , vehicle  velocity. 
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roadway  grade,  and  wind  velocity)  at  each  time  increment.  At  the 
end  of  the  traveling  mode-  (a  known  traveling  time),  the  simulation 
switches  baek  to  Ihe  charging  loop  for  flywheel  recharging.  These 
alternating  loops  are  both  repeal  (id  until  the  desired  number  of  cycles 
of  the  extended  mission  have  been  completed. 

The  program  tabular  outputs  are  prepared  at  each  time  increment  and 
stored  in  the  two  output  data  files.  Listings  of  the  dynamics  and  pow- 
er loss  outputs  during  the  wayside  charging  and  traveling  modes  are 
given  in  Tables  1 and  2,  and  of  the  thermal  outputs  during  both  modes, 
in  Table  3.  Outputs  are  available  in  either  one-second  or  five-second 
time  increments.  Sample  outputs  are  provided  in  Appendix  D,  which 
includes  comment  statements  as  wel]  as  output  data. 

At  the  end  of  each  complete  computer  run,  subroutine  FORM  3 prints 
for  ready  reference  a tabular  summary  of  key  parameters  and  their 
peak  values.  This  summary  appears  at  the  end  of  each  of  the  two  out- 
puts. The  tabular  summary  shown  in  Table  4 is  for  a run  of  five  con- 
secutive mission  cycles,  to  insure  that  the  electric  machine  tempera- 
tures have  stabilized. 

Finally,  the  simulation  provides  a plotting  program.  Chart  F.  4,  that 
displays  significant  system  parameters  in  a 6-channel  .~d  r i p-ehart 
fashion.  A table  (see  Table  3)  describing  the  key  system  parameters 


Table  1 


Dynamic  and  Power  Loss  Outputs  During 


Wayside  Charging  Mode  (FOR  $ 3. 

DAT)* 

Parameter 

Symbol 

Units 

Flywheel  speed 

CO 

o 

pu 

Flywheel  acceleration 

GO 

f 

rad/ 

Alternator  power  input 

p. 

in 

w 

kW 

Alternator  terminal  voltage 

h 

V 

Alternator  armature  current  (per  phase) 

r4 

pu 

Alternator  field  current 

:5 

pu 

Alternator  efficiency 

V 

pu 

Alternator  armature  copper  loss 

A 

1 

kW 

Alternator  field  copper  loss 

A2 

kW 

Alternator  core  loss 

A3 

kW 

Alternator  friction  and  windage  loss 

A4 

kW 

Alternator  stray  load  loss 

A5 

kW 

Alternator  total  loss 

P7 

kW 

Alternator  power  angle 

6 

deg 

Alternator  power  factor  angle 

9 

deg 

* The  order  of  the  output  quantities  in  the  table  corresponds  to  the 
sequence  in  which  they  appear  in  the  computer  tabular  outputs. 
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Table  2 


Dynamic  and  Power  Loss  Outputs  During 

Traveling  Mode  (FOR  </>  3. 

DAT)* 

Parameter 

Symbol 

Units 

Vehicle  velocity 

V 

mi/  h 

Propulsion  (tractive)  force 

F 

lb 

Flywheel  speed 

CO 

o 

pu 

Power  required  from  flywheel 

p 

kW 

Propulsion  power  (at  wheels) 

P2 

kW 

Alternator  power  factor  angle 

9 

deg 

DC  traction  motor  total  loss 

P5 

kW 

Alternator  total  loss 

P7 

kW 

DC  traction  motor  terminal  voltage 

E 

o 

pu 

DC  traction  motor  armature  current 

I 

o 

pu 

Alternator  power  angle 

5 

deg 

Alternator  field  excitation  voltage 

E5 

pu 

Alternator  field  excitation  current 

h 

pu 

Alternator  efficiency 

^ac 

pu 

DC  traction  motor  efficiency 

^dc 

pu 

Solid-state  controller  total  loss 

P6 

kW 

Solid-state  controller  firing  angle 

a 

deg 

Net  propulsion  energy 

E2  to  4 

kWh 

Net  energy  for  all  losses 

E3  to  4 

kWh 

Alternator  armature  current  (per  phase) 

*4 

pu 

* The  order  of  the  output  quantities  in  the  table  corresponds  to 

the  sequence  in  which  they  appear  in  the  computer  tabular  outputs. 


Table  3 


Thermal  Outputs  During 


Wayside  Charging  and  Traveling  Modes 

(FOR  26.  DAT)* 

Parameter 

Symbol 

Units 

DC  traction  motor  total  loss 

P5 

kW 

DC  traction  motor  rotor  input  power 

P 

rdc 

kW 

DC  traction  motor  stator  input  power 

P 

Sdc 

kW 

DC  traction  motor  gap  cooling  power 

pdc 

out 

kW 

DC  traction  motor  rotor  temperature  rise 

T 

rdc 

°C 

DC  traction  motor  stator  temperature  rise 

T 

Sdc 

°c 

DC  traction  motor  gap  temperature  rise 

T 

gdc 

°c 

Solid-state  controller  (PDR)  heating  rate 

P6 

kW 

Solid-state  controller  (PDR)  cooling  rate 

P6 

out 

kW 

Solid-state  controller  (PDR)  junction  temp. 

rise  T, 

D 

t 

P7 

°C 

Alternator  total  loss 

kW 

Alternator  rotor  input  power 

P 

r 

ac 

kW 

Alternator  stator  input  power 

P 

s 

ac 

kW 

Alternator  gap  cooling  power 

P 

ac 

out 

kW 

Alternator  rotor  temperature  rise 

T 

r 

ac 

°C 

Alternator  stator  temperature  rise 

T 

s 

ac 

°c 

Alternator  gap  temperature  rise 

T 

gac 

°c 

* The  order  of  the  output  quantities  in  the  table  corresponds  to  the 
sequence  in  which  they  appear  in  the  computer  tabular  outputs. 
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Table  4 

Tabular  Summary  of  Key  Parameters  and 
Their  Peak  Values* 


. ....  FLYWHEEL  .VEHICLE,  simulation  

V_tRi.CLEL_EO.U_I_V_AL^I_WE:.l  GfcU_s 13  22Z.LRS 

VEHICLE  FRONTAL  AREA  = 53.20  SQ-FT 

ll*  ___45. , 0_5UIG-F  T«S0P  

RESISTANCE  COEFFICIENTS  j 

A = .023003  __  ___  

B = .000175 

C _=_.  *AiO0 0 0 

TOTAL  TRAVELING  TIME  = 1,1319  HRS 

. .TOTAL  DIST  ...TRAVELED  = 20  . 67  MJ. 

PEAK  HEADWIND  = 30,0  MPH 

. PE AK_  V E L 0 C m_= . 3 0_,  0_ftPfci 

PEAK  ACCEL.  = 1,00  MPH/SEC 

PEAK  DECEL.  = -1^00  MPH/SEC 

PEAK  GRADE  = -10*0  PERCENT 

PcAK_IH  RUS_I_= 1903.. UTS 

PEAK  PROPULSION  POWER  = 84.96  KW 

PEAK_DC_L.OSS_s_ ?2 .0?  KW 

PEAK  AC  LOSS  ? 74,04  KW 

PEAK  TOTaL~POHER  a 196.4  KW 

PEAK  DC  ARM,  CURRENT  = 2,186  PU 

PEAK  A C..A.R  M «_C  U R_R  £ NT  = 3,563  PU 

PEAK  AC  POWER  FACTOR  ANGLE  = 180.0  DEG 

PEAK  AC  POWER  ANGLE  = -73.2  DEG 

PEAK  AC  FIELO  CURRENT  = 8,967  PU 

PEAK  A.C^_IXU3— VJOLIAG E_=_ JZ.^8.a5_PU 

PEAK  DC  ROTOR  TEMP  ABOVE  APB.  = 51.3  DEG  C 

...  PEAK  DC  STATOR  TEMP_ABO  VEL_A.MR = 23.6  ..DEG._C_ 

PEAK  AC  ROTOR  TEMP  ABOVE  AHB.  = 321*5  DEG  C 

PEAK. AC  STATOR  TEMP  ABOVE  AM9.  = 372*8  DEG  C 

PEAK  PDR  JUNCTION  T£MP  ABOVE  A MB . = 155.8  DEG  C 

PEAK  HEAT  SINK  TEMP  ABOVE  AMR  = 3.6^.5_DE£_Q_ 


*Total  traveling  time  and  distance  traveled  for  a computer 
run  of  five  consecutive  mission  cycles. 


Table  5 

Key  System  Parameter  Listing 


3-10 


(Masthead  for  Output  Plots) 

FLYUHEEL  VEHICLE  SIMULATION 
PADS* CONFIGURATION  NO.  1 
THERMAL  PLOTS 


VEHICLE  EQUIV.  WT  . 
FRONTAL  AREA 
DRAG  COEFF. 

TIRE  FRICT.  COEFF. 

WIND  VELOCITY 

FLYWHEEL 


R0A0WAY/VEL0CITY 
PROFILE  NO.  1 


TRACTION  MOTOR 


PDR 

ALTERNATOR 


13222.  LB 
53-0  SQ-FT 
0 .610 

ROLLING  = 0.0230 
COULOMB  = 0.000175  PER  MPH 
30.  MPH  (RETARDING)-OUTBOUND  LEG 
-30.  MPH  (AIDING)  - RETURN  LEG 
I = 45.0  SLUG-FT-SQD 


MAX  SPEED 
INITIAL  SPEED 
MAX  SPEED 
TOTAL  DISTANCE 
MAX  GRADE 


12000.  RPM 
6000.  RPM 
30.  MPH 
4 .13  MILES 
10.0  PERCENT 


TOTAL  TRAVEL  TIME  = 0.2264  HR 

70-2  HP  RATING  AT  2730  RPM 

TERM.  VOLT.  = 420.  V 

BASE  ARM.  CUR.  = 140.  A 

FULL  WAVE  3-PHASE  TYPE 

75.8  KVA  AT  12000.  RPM.  3-PHASE 

L-TO-N  TERM.  VOLT-=  215.  V (CONST-) 

BASE  ARM.  CUR-  = 73.2  A/PHASE 


WAYSIDE  RECHARGE  AT  CONSTANT  ARMATURE  CURRENT 


BRAKING 
-0.1018 
-0.4916 
-0.5934 

DATE  OF  RUN  : 28-JUNE-76 


PROPULSION  ENERGY  MOTORING 

WALNUT  TO  ENGR.  STA:  KW-HRS  =•  3.4739 
ENGR.  TO  WALNUT  STA;  KW-HRS  = 1.8396 

ROUND  TRIP;  KW-HRS  = 5.3136 


Propulsion  and  Distribution  System 
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for  a plot  is  also  provided  on  left  side  of  each  output  plot.  The  actual 
plotting  is  done  off-line  by  a Calcomp  plotter.  The  three  following 
sets  of  plots  are  generated,  all  using  mission  time  as  the  abscissa: 

A power  plot  shows:  propulsion  power  during  the 
traveling  mode  and  charging  power  during  the  way- 
side  charging  mode;  power  from  the  flywheel;  and 
total  power  losses  in  the  major  components. 

A dynamic  plot  shows:  the  flywheel  energy;  the  dc 
traction  motor  armature  current;  and  the  solid-state 
controller  firing  angle. 

A thermal  plot  shows  the  temperature  rise  of  the 
hottest  element  of  each  thermally-modeled  component: 
the  armature  of  the  dc  traction  motor;  the  stator  of 
the  alternator;  and  one  thyristor  junction  of  the  solid- 
state  controller. 

To  facilitate  comparisons  among  these  three  plots,  the  first  three 
channels  of  each  plot  show  the  vehicle  velocity  profile,  roadway  grade 
profile,  and  computed  tractive  force.  Samples  of  these  plots  are 


shown  in  Figs.  18  thru  21. 


If  flywheel-spin-down,  which  results  either  from  an  undersized 
flywheel,  or  from  too  much  drain,  does  not  permit  the  alternator 
to  generate  the  required  ac  voltage,  a special  feature  of  the  simu- 
lation stops  the  program  and  prints  an  error  message  that  prevents 
further  useless  computation. 

Calculation  of  Propulsion  Power  and  Auxiliary  Power  Requirements 

The  total  power  required  from  the  flywheel  at  any  instant  of  travel 
time  is  the  sum  of  three  requirements: 

Propulsion  power 
Auxiliary  power 
Total  power  losses 

Because  the  detailed  simulation  of  the  loss  process  for  all  of  the 
major  components  represents  a basic  departure  from  a conventional 
approach,  it  will  be  described  in  a separate  section.  The  methods 
for  computing  the  auxiliary  and  propulsion  power  requirements  at 
each  instant  of  time  are  more  conventional.  Since  auxiliary  power 
is  the  least  significant  power  requirement,  it  warrants  only  brief 


discussion. 


The  required  auxiliary  power  has  been  estimated  in  the  simulation. 
Its  value  can  vary  widely  as  a function  of  the  ambient  temperature 
of  the  day  and  related  heating  or  cooling  requirements  of  the  vehicle 
interior.  The  requirement  is  modeled  as  a constant  3 -phase  load  of 
6 kW,  0.  8 PF  lag,  during  the  traveling  and  charging  modes.  It  is 
supplied  from  the  alternator  during  motoring,  from  the  solid-state 
controller  during  regenerative  braking,  and  from  the  wayside  power 
connector  during  wayside  charging. 

The  methods  used  to  determine  the  required  propulsion  power  are 
well-known.  The  propulsion  power  at  any  instant  of  travel  time  is 
the  product  of  the  tractive  force  and  the  /ehicle  velocity.  Velocity 
is  a known  input  to  the  simulation.  Tractive  force  is  calculated  as 
the  sum  of  four  component  forces: 

Inertia 
Gravity 
Tire  friction 
Aerodynamic  drag  , 

Inertia  forces  are  the  forces  required  to  change  vehicle  velocity  to 
conform  to  the  velocity  profile.  They  are  the  product  of  the  vehicle 
equivalent  mass  and  the  slope  of  the  vehicle  velocity/time  profile 


(acceleration).  The  vehicle  equivalent  mass  is  the  sum  of  the  actual 


vehicle  mass  and  the  referred  value  of  the  rotary  inertia  of  those  ro- 
tating parts  whose  speed  is  directly  proportional  to  vehicle  speed, 
such  as  the  drive-wheels  and  the  dc  traction  motor.  For  example, 
for  an  equivalent  vehicle  mass  of  13,  222  lb  and  an  acceleration  of 
1.4  mi/h/s,  the  inertia  force  is  884  lb. 

Gravity  forces  equal  the  product  of  the  vehicle  equivalent  weight  and 
the  sine  of  the  roadway  grade  angle  above  the  horizontal.  For  the 
same  vehicle  on  a 5 percent  grade,  the  gravity  force  is  661  lb. 

Tire  friction  forces  result  from  the  resistance  of  the  roadway  to  tire 
motion.  They  depend  on  vehicle  weight  and  velocity,  tire  pressure, 
and  tire  design.  The  tire  friction  model  used  here,  which  is  based 
on  nominal  conditions,  is: 


F = uW  (1  + C V)  (MF)  lb 
tire  c 


tire 


where: 


rolling  friction  coefficient  (0.0230  lb/ lb) 


W 


C 


V 


c 


vehicle  weight  (12,  238  lb) 

_4 

coulomb  friction  coefficient  (1.  75  x 10  h / mi) 

vehicle  velocity  (ft/ s) 


MF  = mass  factor-ratio  of  equivalent  to  actual  mass  . 
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For  example,  for  a vehicle  velocity  of  30  mi/h,  the  tire  friction 
force  is  306  lb. 


Aerodynamic  drag  forces  result  from  the  air  resistance  to  vehicle 
velocity  and  wind  velocity.  For  this  simulation,  the  wind  speed  is 
taken  as  constant  (30  mi/h)  in  a vehicle  frame  of  reference.  The 
wind  speed  opposes  vehicle  motion  during  the  uphill  outbound  seg- 
ment of  the  trip  and  aids  it  during  the  downhill  inbound  segment.  The 
drag  force  model  is: 


drag 


1/2  p (V+  V ) 
w 


V + V 


w, 


V + V I C,A 


lb 


w 


where : 


p = air  density  at  sea  level  (0.  002378  slugs/ ft  ) 

V = vehicle  velocity  (ft/  s) 

V = wind  velocity  (ft/  s)  (negative  for  tailwind) 
w 

C = vehicle  drag  coefficent  (0.  6l) 
d 

2 

A = vehicle  frontal  area  (53.  0 ft  ). 


For  example,  for  the  same  vehicle  at  a velocity  of  30  mi/h  and  a head- 
wind of  30  mi/h,  the  aerodynamic  drag  force  is  298  lb. 


The  sum  of  the  force  components  for  the  13,  222  lb  vehicle  accelerating 
at  1. 4 mi/h/s,  at  a velocity  of  30  mi/h,  ona  5 percent  grade  and  against  a 


headwind  of  30  mi/h,  is  2,  109  lb.  A flowchart  for  computing  the 
tractive  force  and  propulsion  power  is  shown  in  Fig.  9. 
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*If  grade  profile  input  is  provided  in  the  form  of  G (x)  rather  than  G (t), 
then  G (t)  = G (x)^  x V (t)„ 


Fig.  9 

Flowchart  for  Computing 
Tractive  Force  F and  Propulsion  Power  P 
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4,  MODELS 


Procedure 

Calculation  of  the  power  losses  in  the  energy-storage  unit  and  the 
power  train  is  required  in  the  simulation  in  order  to  compute  the  power 
delivered  by  the  flywheel  and  the  temperatures  of  the  components 
at  each  instant  of  travel  time.  The  power  loss  for  each  component 
is  calculated  by  a two-step  process,  using  first  the  dynamic  model, 
then  the  power  loss  model  for  the  component.  The  power  loss  model 
utilizes  as  inputs  the  electrical  variables  and  speeds  which  are  the 
outputs  of  the  dynamic  model.  The  temperature  rises  for  the  elec- 
tric machines  and  the  thyristors  in  the  controller  are  calculated  in  a 
third  step,  using  thermal  models  for  the  components  and  the  losses 
from  the  power  loss  models.  These  models  are  derived  in  this 
section;  numerical  values  for  the  parameters  are  given  for  the  com- 
ponents of  the  baseline  system.  In  order  to  use  the  simulation  as  a 
design  tool,  and  to  vary  the  machine  ratings,  the  relationships  are 
given  in  terms  of  the  linear  dimension  of  the  electric  machines  between 
the  power  ratings,  loss  components  and  electrical  parameters. 
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Dynamic  and  Power  Loss  Models 

The  dynamic  models  of  each  component  of  the  electric  propulsion 
system  are  basically  the  equations  that  relate  the  input  and  output 
variables.  The  dynamic  models  of  the  electric  machines  are  more 
complicated  than  the  other  components.  The  models  are  sets  of 
equations  that  relate  the  speed  and  torque  of  the  electric  machines 
to  the  electrical  variables  of  current  and  voltage.  These  equations 
are  based  on  the  machine  parameters  of  resistance,  reactance,  and 
saturation,  which  are  either  obtained  from  the  manufacturer  or  cal- 
culated from  the  per-unit  quantities  for  machines  of  the  same  type 
and  size. 


The  power  loss  models  are  sets  of  equations  that  relate  the  elec- 
trical variables  and  speeds  of  the  major  components  to  their  power 
losses.  These  models  are  generally  more  elaborate  for  the  baseline 
system  components  that  are  expected  to  exhibit  the  largest  losses, 
i.  e.  , the  lowest  efficiencies.  The  most  elaborate  loss  models,  each 
including  five  separate  loss  terms,  are  for  the  alternator  and  the  dc 
traction  motor.  A two-term  loss  model  describes  the  solid-state 
controller  losses.  The  model  that  represents  the  transmission  losses  is 
simplest,  since  a constant  percentage  of  the  input  is  assumed  lost  in  that 
component  power.  The  flywheel  mechanical  losses  (friction  and  windage) 
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were  not  modeled,  since  test  data  for  representative  flywheels 
running  in  evacuated  enclosures  were  not  available,  but  could  be 
added  at  a later  date  if  necessary.  Flywheel  losses  are  expected 
to  be  small  and  could  be  lumped  with  the  alternator  friction  and 
windage  losses. 

As  the  flywheel  is  spun-up  in  the  wayside  charging  mode,  the  only 
losses  occur  in  the  alternator,  which  is  used  as  a motor,  while  the 
solid-state  controller,  dc  traction  motor,  and  transmission  are  not 
in  operation.  The  auxiliary  load  is  assumed  to  be  supplied  from 
the  wayside  power  connector.  During  the  traveling  mode,  regard- 
less of  the  direction  of  power,  all  four  major  components  develop 
losses. 

Part  of  the  output  of  the  simulation  is  a tabulation  of  all  of  the 
electrical  variables,  machine  speeds,  loss  terms,  component  losses, 
and  total  losses.  In  a further  calculation,  the  component  losses  are 
used  with  the  thermal  models  to  obtain  the  temperatures  of  the  com- 
ponents at  each  instant  of  travel  time. 

The  dynamic  and  loss  models  are  developed  in  this  section  in  the 


order  of  calculation  at  each  time  increment  of  the  computer  program, 
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that  is,  from  the  drive  wheels  through  the  propulsion  system  to  the 
flywheel.  Hence  the  dc  traction  motor  model,  which  effectively 
includes  the  gear  box  is  addressed  first,  followed  by  the  controller 
and  the  alternator.  The  alternator  is  the  most  complicated  of  the 
component  models. 

DC  Traction  Motor.  The  equations  for  the  dynamic  model  of  the 
dc  traction  motor  are  given  in  Appendix  A.  Basically,  the  armature 
current  is  computed  from  the  tractive  force  at  the  drive-wheels  as 
reflected  through  the  gearbox.  The  motor  speed  is  also  related  to 
the  drive-wheel  speed. 


The  equations  can  be  written  in  per-unit  terms  or  in  terms  of  the 
variables  in  ohms,  volts,  amps,  r/  min,  etc.  The  per-unit  terms 
are  more  convenient  because  the  value  of  the  variable  conveys  its 
proportion  to  the  rating  of  the  machine;  1^  = 0.  5 pu,  for  example, 
means  that  the  armature  current  is  at  half  its  base  value.  An 
additional  advantage  is  that  for  small  changes  of  the  power  rating  of 
the  machine,  the  per-unit  parameters  hardly  change;  the  same  equa- 


tions can  usually  be  used. 
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Before  the  per-unit  equations  for  the  dynamic  model  can  be  written, 
the  base  quantities,  corresponding  to  1. 0 pu,  must  be  given  for  the 
machine.  Usually,  these  quantities  are  the  nameplate  rated  quanti- 
ties. The  base  quantities  for  the  dc  traction  motor  used  in  the  base- 
line system  are  given  in  Table  6. 


Table  6 


DC  Traction  Motor  Base  Quantities 


Armature  current 

I 

o 

= 140  A 

Armature  voltage 

Vdc 

= 420  V 

Armature  resistance 

R 

adc 

= 3.0  <T2 

Speed 

N 

= 2730  r/  min 

Torque 

T 

= 135  lb/ ft 

Mechanical  power 

P 

m 

= 52.  3 kW  (70  hp) 
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The  five  equations  for  the  power  loss  model  are  written  in  terms  of 
the  per-unit  variables  but  yield  the  power  loss  components  in  kW. 
The  equations  for  the  armature  circuit  copper  and  independent  field 
losses  (0^  and  D^)  are  derived  from  the  circuit  resistances.  The 
equation  for  friction,  windage,  and  iron  losses  (D^)  assumes  a 
variation  with  speed  to  the  2.  5 power.  The  stray  load  loss  (D^)  is 
taken  as  1 percent  of  the  mechanical  power.  The  brush  loss  (D  ) is 
derived  from  a constant  2-V  brush  drop.  The  five  equations  in  per- 
unit  variables  yield  the  following  loss  components  in  kW : 


D! 

= 4.  10  I 

a 

- Armature  circuit  copper  loss 

Dz 

- 0.  288  (N)2-  5 

- Independent  field  loss 

°3 

= 1.  56  (N/  2730) 

- Friction,  windage,  and  iron  losses 

°4 

= io~2  P * 
m 

- Stray  load  loss 

°5 

ii 

o 

tv 

oc 

o 

o 

- Brush  loss  , 

In  the  computer  program,  the  variable  Pm(in  D ) is  expressed 
as  P /n  for  motoring  and  P x q for  braking,  where  P? 


2'  "gb  *”*  & "2  gb 

is  the  propulsion  power  at  the  wheels  and  q 

efficiency  (0.  92) . 


gb 


is  the  gearbox 
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The  sum  of  the  loss  components  calculated  at  the  base  conditions 
from  the  power-loss  equations  must  yield  the  losses  at  the  base  con- 
ditions calculated  from  the  efficiency.  The  loss  components  for  the 
baseline  system's  dc  traction  motor  are  shown  in  Table  7.  The 
manufacturer  usually  supplies  the  armature  resistance,  field  resis- 
tance, and  full-load  efficiency. 


Table  7 


DC  Traction  Motor  Loss  Components  at  Base  Conditions 


D, 

= 4.  10  kW 

- Armature  circuit  copper  loss 

°2 

= 0.  29  kW 

- Independent  field  loss 

°3 

= 1.  56  kW 

- Friction,  windage,  and  iron  losses 

°4 

- 0.  52  kW 

- Stray  load  loss 

°5 

= 0.  28  kW 

- Brush  loss 

P5 

= 6.  75  kW 

- Total  loss 

The  stray  load  and  brush  loss  components  can  be  found  as  previously 
described.  The  friction,  windage,  and  iron  loss  component  makes  up 
the  remainder  of  the  losses,  which  can  be  calculated  from  the  efficiency, 
i.  e.  once  the  efficiency  is  known. 
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Solid-State  Controller.  The  equations  for  the  dynamic  model  of  the 
controller  relate  the  dc  output  voltage  and  current  to  the  ac  input 
voltage  and  current  and  yield  the  controller  firing  angle  as  a depen- 
dent variable.  The  power  loss  model  does  not  require  the  results 
of  the  dynamic  model  calculation.  The  power  losses  consist  of  two 
components:  the  thyristor  forward  drop  losses  and  the  snubber  network 
losses.  Forward  drop  loss  of  this  6-thyristor  controller  are  based  on 
typical  values.  Forward  drop  loss  is  assumed  to  vary  with  direct 
current  and  is  taken  at  0.  3 kW/  100  A,  i.  e. , 3 V drop.  Snubber  network 
losses  are  assumed  to  be  1 percent  of  the  output  power  of  the  solid-state 
controller.  When  the  dc  traction  motor  delivers  its  rated  output,  the 
sum  of  both  loss  terms  is  roughly  1 kW. 

« 

Alternator.  The  dynamic  model  for  the  alternator  is  much  more 
complicated  than  the  model  for  the  traction  motor,  for  the  following 
reasons.  First,  since  the  alternator  is  an  ac  machine,  the  stator 
electrical  quantities  must  be  expressed  in  complex  numbers  to 
identify  magnitude  and  phase  of  the  currents  and  voltages.  Second, 
the  baseline  alternator  operates  over  a range  of  speeds,  typically 
from  0.  5 to  1.0  pu;  practically  all  conventional  alternator  analysis  is 
based  on  constant-speed  operation.  Third,  to  obtain  a sufficiently 
accurate  dynamic  model,  the  alternator  must  be  represented  in  its 
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direct  and  quadrature  (d-q)  axes.  Fourth,  because  the  alternator 
will  operate  down  to  a speed  of  0.  5 pu  and  at  high  field  current  to 
maintain  terminal  voltage,  saturation  of  the  magnetic  circuit  must 
be  included  in  the  model. 

The  equations  for  the  dynamic  model,  which  are  given  in  Appendix  A, 

use  the  standard  machine  impedance,  d-q  axis  representation,  and 

11,12 

the  Kingsley  method  for  saturation  correction.  They  describe 

the  operation  of  the  alternator  in  the  traveling  mode,  as  the  flywheel 
and  alternator  speed  declines  from  1.  0 pu  to  as  low  as  0.  5 pu;  or 
in  the  wayside  charging  mode,  with  the  alternator  operating  as  a 
synchronous  motor  as  the  speed  rises  from  0.  5 pu  to  1.0  pu.  In  the 
traveling  mode,  the  alternator  delivers  energy  to  the  solid-state 
controller  and  supplies  the  auxiliary  load;  the  terminal  voltage  is 
held  constant  over  the  speed  range.  In  the  wayside  charging  mode, 
the  alternator  receives  power  from  the  wayside  inverter. 

The  per-unit  system  used  for  the  dynamic  model  defines  the  currents, 
voltages,  and  reactances  at  the  base  ojq  = 1.0  pu.  At  other  than  base 
speed,  the  per-unit  reactance  is  multiplied  by  the  per-unit  speed  to 


reflect  the  effect  of  the  frequency  on  the  reactances.  If  the  per-unit 
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impedances  are  not  provided  by  the  manufacturer  for  a specific 

alternator,  typical  per-unit  values  can  be  obtained  for  alternators 

1 3 

of  the  same  similar  rating  and  speed  from  a handbook.  The  per-unit 
impedances  for  the  baseline  alternator  --  a 75-kVA,  1 2,  000-r/  min, 
4-pole  unit  --  have  been  provided  by  Westinghouse  and  are  given  in 
Table  8. 


fable  8 


Alternator  Per-Unit  Impedances 


r 

a 

ac 

- 0.  0495 

- Resistance  (200°C) 

Xd 

= 3.  40 

- D-axis  unsaturated  reactance 

X 

= 1. 24 

- Q-axis  reactance 

q 

X 

a 

= 0.085 

- Leakage  reactance 

1 1 

= 0.  145 

- Subtransient  reactance 

d 


Ratings  - The  power  and  speed  rating  of  the  alternator  determine 
its  physical  size  and  its  ability  to  carry  continuous  load.  Selection  of  an 
alternator  should  be  based  on  three  criteria:  (1)  its  ability  to  provide 
peak  power  at  some  point  in  the  mission;  (2)  its  energy  losses  between 
wayside  charging  stations;  and  (3)  its  maximum  internal  temperature 


as  compared  to  the  allowable  temperature.  For  an  alternator 
running  at  one  speed  and  variable  load,  the  rating  can  be  selected 
by  calculating  the  rms  current  over  an  operating  period.  However, 
for  an  alternator  operating  over  a range  of  speed  and  load,  as  in 
a flywheel  energy-storage  vehicle,  the  selection  can  best  be  made  by 
using  the  dynamic,  loss,  and  thermal  models  described  in  this  report 

After  power  and  speed  ratings  are  determined,  the  alternator  voltage 
rating  must  be  selected.  For  a given  power  rating  and  frame,  the 
voltage  and  current  ratings  of  the  machine  are  like  scale  factors. 

Any  values  can  be  used,  provided  their  product  yields  the  given 
power  rating.  However,  the  ratings  selected  will  significantly  affect 
the  alternator  power  losses.  In  the  traveling  mode,  the  alternator  is 
assumed  to  operate  at  constant  terminal  voltage,  while  the  speed 
ranges  from  1. 0 pu  to  0.  5 pu.  If,  on  one  hand,  the  rated  voltage  is 
set  equal  to  the  required  voltage,  then  at  0.  5 pu  speed  the  alternator 
must  be  capable  of  twice  the  air-gap  flux  density  developed  at  1 . 0 pu 
speed.  The  field  current  and  field  losses  will  be  high  at  low  speed. 

If,  on  the  other  hand,  the  rated  voltage  is  set  equal  to  twice  the 
required  voltage,  then  the  increase  of  stator  turns  will  raise  the  arma 
ture  resistance  and  make  the  armature  losses  high  at  all  speeds.  In 
the  baseline  alternator,  the  rated  voltage  was  expressed  as  a rewind 
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factor  (RWF)  times  the  original  115  V/ phase  rating  of  the  machine. 
A rewind  factor  of  3.  0 was  selected. 

The  base  quantities  for  the  alternator  used  in  the  baseline  system 
are  shown  in  Table  9. 


Table  9 

Alternator  Base  Quantities 


Power 

P 

= 75.  8 kVA 

Voltage 

vt 

= 345  V (i-n) 

Armature  current 

h 

= 73. 2A 

Speed 

u 

o 

= 12,  000  r/min 

Field  current 

h 

= 14.  3A 

Base  impedance 

z. 

o 

= 4.  70ft 

The  saturation  curves  for  the  baseline  alternator  are  shown  in  Fig.  10 
for  w = 1.  0 pu  and  0.  5 pu.  The  saturation  curve  for  uq  = 1.0  pu  was 
supplied  by  the  manufacturer.  The  saturation  curve  for  ioq  0.  5 pu, 
or  for  any  other  speed,  is  merely  scaled  down  vertically  by  pu  for 
each  value'  of  field  current.  For  the  alternator  under  load,  the  satura- 
tion curve  is  assumed  to  represent  the  air-gap  voltage  vs.  the  net 
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*Air  gap  voltages  (E  ) for  dc  traction  motor  operating  at  rated 
power  and  speed. 


Fig.  10 


Saturation  Curves  for  Alternator 


L 
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direct-axis  ampere  turns  expressed  in  equivalent  field  amperes.  At 
no-load  condition,  of  course,  the  air-gap  voltage  is  identical  with  the 
terminal  voltage;  the  net  direct-axis  ampere  turns  are  merely  the  field 
ampere  turns. 

The  air-gap  voltage  point  represents  the  degree  of  saturation  of  the 

machine.  For  base  operation  with  PF  =1.0,  the  air-gap  voltage 

is  E = 362  V (Point  O).  For  conditions  when  the  dc  traction  motor 
r 

is  operating  at  rated  70  hp  and  rated  2730  r/  min  while  the  alternator 
is  running  at  w =1.0  pu,  the  air-gap  voltage  is  = 265  V (Point  M); 
i.  e.  , in  other  words,  the  alternator  is  less  saturated  than  it  is  during 
base  operation.  For  the  same  load,  but  with  the  alternator  running  at 
= 0.  5 pu,  the  alternator  is  heavily  saturated  at  E^=  252  V (Point  N). 

Phasor  diagrams  - The  dynamic  model  of  the  alternator  is 
expressed  as  a phasor  diagram.  It  is  used  to  calculate  air-gap  vol- 
tage E^  and  the  field  current  1^  at  each  instant  of  time  during  the 
mission  to  obtain  two  of  the  five  loss  components.  The  phasor  dia- 
gram for  the  alternator  at  base  conditions  is  shown  in  Fig.  11.  The 
air-gap  voltage  E^  is  calculated  first  to  find  the  degree  of  saturation; 
the  position  of  the  direct  axis,  identified  with  the  power  angle  6 , is 

calculated  next  using  the  quadrature -axis  reactance  x . The  magne - 

q 

tizing  portion  of  the  direct-axis  reactance  (x  -xj  is  corrected  for 
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0 


2.33  pu 


1.24  pu 


= 0.085 
0.05  pu 


pu 


* Alternator  base  conditions:  345  V (^-n);  75.  8 kW;  1.  0 PF;  12,  000 

Fig.  11 

Alternator  Phasor  Diagram,  Traveling  Mode: 
Generator  Operation  at  Alternator  Base  Conditions* 


/min. 
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saturation  using  the  slope  of  the  Kingsley  air-gap  line  through  the 
air-gap  voltage  point,  as  shown  in  Fig.  10.  The  diagram  is  finally 
completed  to  find  the  excitation  voltage  E To  find  the  field  current 
in  amperes,  the  per-unit  excitation  voltage  must  be  corrected  for 
saturation  by  using  the  ratio  of  the  slopes  of  the  Kingsley  air-gap 
line  relative  to  the  base  air-gap  line.  For  example: 

Corrected  E = E / slope  ratio  = 2.  33  pu / 0.  954  = 2.  45  pu 
5 5 

I = 2.  45  pu  x 14.  3 A = 35  A . 

5 

The  phasor  diagrams  for  the  alternator  at  two  points  in  the  traveling 
mode  are  shown  in  Fig.  12.  These  diagrams  correspond  to ‘the  con- 
ditions of  the  flywheel  and  alternator  at  full  speed  and  half  speed 
while  the  dc  traction  motor  is  running  at  its  base  conditions.  The 
operating  points  are  also  shown  on  the  saturation  curves  of  Fig.  10. 

The  terminal  voltage  is  0.  62  pu  or  215  V;  the  armature  current  is 
1 . 65  pu  or  1 2 1 A,  at  PF  = 0.  85  lag. 

As  Fig.  10  shows,  at  full  speed,  point  M,  the  alternator  is  hardly 
saturated;  the  field  current  is  3.  91  pu,  or  56  A.  At  half  speed,  point 
N,  the  alternator  is  heavily  saturated.  Since  the  field  current  is  9.  95  pu, 
or  142  A,  both  the  armature  and  field  losses  are  extremely  high.  Over 
a typical  mission,  the  alternator  might  have  to  carry  even  higher 
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jl,  X,  W = 5 . 2 2pu 
4f  dsat  0 


E r=  5 . 5 3pu 


6=  4 
0 = 3 


Vt, 


Vt, 


3.  1 6 


pu 


u.  49  pu 


sat. 


sa 


(a)  Alternator  at  Full  Speed 


(b)  Alternator  at  Half  Speed 


^Alternator  Operation:  215  V (f-n);  66.  1 kW;  0.  85  PF  lag;  12,  000  and  6,  000  r/min. 

Fig.  12 

Alternator  Phasor  Diagrams,  Traveling  Mode: 


Generator  Operation  at  DC  Traction  Motor  Base  Conditions- 


armature  currents  to  support  the  tractive  force  requirements  of  the 

dc  traction  motor.  If  such  high  armature  currents  occur  when  the 

flywheel  has  run  down  to  half  speed,  the  alternator  field  current  might 

be  even  higher  than  the  9.  95  pu  value,  and  increased  losses,  propor- 
2 

tLonal  to  I_  would  result, 
o 

The  phasor  diagrams  for  the  alternator  operating  as  a synchronous 
motor  at  points  in  the  wayside  charging  mode  are  shown  in  Fig.  13. 
These  diagrams  are  based  on  2.  2 pu  armature  current  at  PF  =1.0. 
The  terminal  voltage  increases  linearly  with  speed  to  a final  value 
of  1.2  pu  when  the  alternator  is  at  rated  speed.  During  the  charging 
period,  at  half  speed,  the  alternator  is  not  saturated;  the  field  cur- 
rent is  5.  02  pu,  or  71.  8 A.  At  full  speed,  the  alternator  is  partially 
saturated;  the  field  current  is  5.  18,  or  74.  0 A.  The  field  current  is 
practically  constant  during  this  mode  of  operation.  The  alternator 
operating  as  a motor  could  operate  at  even  higher  power,  and  thus, 

shorten  the  charging  time  if  the  final  value  of  the  terminal  voltage 
were  raised  above  1.2  pu.  However,  the  field  current  and  losses,  as 

well  as  the  temperature  rise,  would  be  higher  than  for  the  assumed 
charging  condition.  The  limitation  on  charging  is  the  pull-out 
power  for  the  alternator,  operating  as  a synchronous  motor. 
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♦Input  to  Machine:  At  half  speed,  207  V,  100  kW;  at  full  speed,  414  V,  200  kW. 

Fig.  13 

Alternator  Phasor  Diagrams,  Wayside  Charging  Mode: 


= 0.  llpu 
.19  pu 


2.73° 
•18  pu 


Alternator  Operating  as  a Synchronous  Motor* 
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Power  loss  model  - The  five  equations  for  the  power  loss  model 
are  written  in  terms  of  the  per-unit  variables  but  yield  the  power 
loss  components  in  kW.  The  equations  for  the  armature  and  field 
circuit  copper  losses  (A^  and  A^)  are  derived  from  the  circuit  resis- 
tances,  including  that  of  the  exciter.  The  equation  for  the  iron  loss 
(A^)  assumes  a variation  with  air-gap  voltage,  to  the  2.  0 power. 

The  air-gap  voltage  is  proportional  to  the  product  of  speed  and 

flux.  The  friction  and  windage  loss  (A^)  is  assumed  to  vary 
3 

as  . The  stray  load  loss  (A^)  has  a constant  term  and  a term 
dependent  upon  armature  current.  The  five  equations  in  per-unit 
variables  for  the  baseline  alternator  yield  the  following  loss  compon- 
ents in  kW: 


A,  =3.94I42 

2 

A_  = 0.  280  + 0.  24  I 

2 5 

A_  = 1. 45  E 2 

3 r 

A . = 2.  05  w 3 

4 o 

An  = 0.  207  + 1.  67  T 2 

5 4 


Armature  copper  loss 

Exciter  and  field  circuit 
copper  loss 

Iron  loss 

Friction  and  windage  loss 
Stray  load  loss  . 


At  the  base  conditions,  the  sum  of  the  loss  components  given  by  the 
power-loss  equations  must  equal  the  known  total  loss  of  the  alternator. 
The  loss  components  for  the  alternator  are  shown  in  Table  10. 
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Table  10 

Alternator  Loss  Components  at  Base  Conditions 


A1 

3.  94  kW 

Armature  copper  loss 

Az 

= 

1 . 63  kW 

- 

Exciter  and  field  circuit  copper  loss 

A3 

= 

1.  60  kW 

- 

Iron  loss 

A4 

= 

2.  04  kW 

- 

Friction  and  windage  loss 

A5 

= 

0.  39  kW 

- 

Stray  load  loss 

P7 

= 

9.  60  kW 

- 

Total  loss 

The  manufacturer  usually  supplies  the  armature  and  field  circuit 
resistances  and  the  full-load  efficiency.  Once  efficiency  is  known, 
the  total  loss  in  kW  can  be  calculated.  The  loss  components  A , A^, 

and  A can  be  allocated  using  typical  values  from  a handbook  or  text 
on  electric  machine  design. 

The  alternator  losses  during  the  traveling  mode,  when  the  dc  traction 
motor  is  operating  at  its  base  condition  (2730  r/ min  and  70  hp) , are 
shown  in  Table  11  for  two  conditions:  when  the  alternator  is  at  full 
speed  and  at  half  speed.  These  operating  points  correspond  to  the 
phasor  diagrams  in  Fig.  12.  Because  the  armature  current  is  1.  65  pu. 
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both  the  armature  and  the  field  losses  are  higher  than  for  the  base 
conditions  of  Table  10.  Compared  to  the  base  loss  of  only  9.  60  kW, 
the  total  loss  varies  from  19.  4 kW  at  full  speed  to  37.  5 kW  at  half 
speed.  The  relative  magnitude  of  the  total  loss  at  full  speed  and  at 
half  speed  can  be  shifted  by  selecting  the  base  voltage  for  the  alternator. 
Selecting  a base  voltage  higher  than  345  V will  reduce  the  field 
loss  A^  but  raise  the  armature  resistance  and  thus  the  armature  loss 
A . Selecting  a base  voltage  lower  than  345  V will  produce  the  oppo- 
site effects  on  the  loss  components  A and  A^. 


Table  11 


Alternator  Loss  Components  During  Traveling  Mode* 


Loss  Component 

Alternator 

Full-speed 

Half-speed 

Ai 

11.  6 kW 

1 1 . 6 kW  - 

Armature  copper  loss 

A2 

4.  0 kW 

24.  0 kW  - 

Exciter  and  field  circuit  copper 

A3 

0.  9 kW 

0.  8 kW  - 

Iron  loss 

A4 

2.  1 kW 

0.  3 kW  - 

Friction  and  windage  loss 

A5 

0.  8 kW 

0.  8 kW  - 

Stray  load  loss 

P7 

19.  4 kW 

37.  5 kW  - 

Total  loss 

>jc 


Constant  Load  of  78.  0 kVA,  66.  1 kW,  215  V ( £ -n) 
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The  alternator  losses  during  wayside  charging  are  shown  in  Table  12 
for  the  alternator  at  half  speed  and  at  full  speed.  Compared  to  the 
base  losses  of  only  9.  60  kW,  the  total  losses  are  practically  con- 
stant, ranging  from  27.  2 kW  to  30.  8 kW. 

Table  12 

Alternator  Losses  During  Wayside  Charging  Mode* 


Loss  Component 

Alternator 

Half-speed 

Full-speed 

Ai 

19.  1 kW 

19.  1 kW  - 

Armature  copper  loss 

A2 

6.  3 kW 

6.  7 kW  - 

Exciter  and  field  circuit  copper 

A3 

0.  4 kW 

1.  8 kW  - 

Iron  loss 

A4 

0.  3 kW 

2.  1 kW  - 

Friction  and  windage  loss 

A5 

1.  1 kW 

1.  1 kW  - 

Stray  load  loss 

P7 

27.  2 kW 

30.  8 kW  - 

Total  loss 

^Armature  current  2.2  pu;  PF  = 1.  0 (motoring);  terminal  voltage 
proportion  to  speed,  up  to  1.  2 pu  at  full  speed. 
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Propulsion  System  Losses.  The  combination  of  the  dynamic  and 
power  loss  models  for  the  energy-storage  unit  and  power-train  com- 
ponents will  yield  both  the  component  and  the  overall  propulsion 
system  power  losses  at  each  instant  of  time  during  the  mission  cycle. 
For  example,  the  component  losses  when  the  dc  traction  motor  is 
running  at  its  base  condition  and  the  flywheel  is  at  full  speed  are 
shown  in  Fig.  14.  The  dc  traction  motor  is  delivering  70  hp,  or 
48  kW,  at  the  drive-wheels.  The  flywheel  must  supply  85.  3 kW  to 
overcome  the  alternator  and  power-train  losses  and  provide  the 
required  auxiliary  load  of  6 kW.  The  overall  propulsion-system 
efficiency  is  62  percent.  If  the  flywheel  were  running  at  less  than  full 
speed  and/  or  the  dc  traction  motor  providing  more  than  its  base  torque 
of  135  lb  ft,  the  propulsion-system  losses  would  be  even  greater  than 


those  under  base  conditions. 
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Thermal  Models 


The  temperatures  of  the  major  components  of  the  drive  system  at 
each  instant  of  travel  time  are  calculated  by  using  thermal  models. 

The  power  losses  required  for  the  calculation  are  obtained  from  the 
power  loss  models.  The  maximum  component  temperatures  show 
whether  the  components  have  been  properly  selected  to  withstand  the 
service  conditions  of  the  vehicle. 

The  insulations  of  rotating  electric  machines,  i.  e.  , the  alternator 
and  dc  traction  motor,  are  assigned  maximum  allowable  hot  spot 
temperatures  by  the  manufacturer.  Whenever  the  operating  temper- 
ature exceeds  the  allowable  temperature,  the  service  life  of  the 
machine  is  reduced,  typically  by  a factor  of  two  for  each  10°C  that 
the  limit  is  exceeded.  If  temperature  excesses  are  large,  e.  g.  , 
more  than  60° C,  complete  failure  of  the  component  may  result.  If 
the  machine  does  not  reach  its  temperature  limit,  even  during  heav- 
iest conditions  of  load,  the  machine  has  been  selected  too  conservative 
and  will  be  too  heavy.  One  purpose  of  the  simulation  is  to  provide  the 
designer  with  means  for  selecting  optimal  ratings  for  the  major  com- 


ponents of  the  drive  system. 
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Unlike  the  rotating  electric  machines,  the  thyristors  in  the  solid-state 
controller  will  burn  out  if  the  junction  temperature  exceeds  the  allow- 
able value  even  for  a fraction  of  a second.  The  temperature  rise  can 
be  reduced  by  selecting  larger  thyristors  with  greater  thermal  mass  and 
conductivity  and/or  selecting  thyristor  heat  sinks  with  greater  thermal 
mass  and  higher  heat  dissipation  rates. 

Thermal  models  were  not  developed  for  the  gearbox  nor  the  flywheel. 
The  gearbox  losses  are  small  and  its  temperature  rise  is  therefore 
not  expected  to  be  large.  Furthermore,  because  the  gearbox  has  a 
large  heat  capacity,  the  effects  of  overheating  would  not  be  drastic. 

The  flywheel  does  not  require  a thermal  model  to  calculate  the  tem- 
perature of  its  bearings.  They  must  be  selected  to  withstand  full 
speed  regardless  of  the  mission  of  the  vehicle. 

The  thermal  models  use  lumped  parameters  of  thermal  mass  and 
thermal  conductivity  to  represent  the  actual  distributed  parameters 
of  the  electric  machines  and  the  thyristors  of  the  solid-state  controller. 
The  power  losses  are  the  independent  variables  (inputs);  the  tempera- 
tures are  the  dependent  variables  (outputs).  Since  the  temperature 
variation  lags  in  time  the  power  loss  inputs,  vehicle  operation  must 
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be  simulated  long  enough  --  i.  e.  , over  several  mission  cycles  - -to 
allow  the  temperature  variations  to  reach  a steady-state  pattern. 

Rotating  Electric  Machines.  Similar  thermal  models  are  used  for 
the  alternator  and  dc  traction  motor.  The  rotor  is  represented  by 
one  thermal  mass  and  the  stator  by  a second  thermal  mass.  Addi- 
tional lumped  elements  can  be  used  at  the  discretion  of  the  modeler; 
two  thermal  masses  are  sufficiently  accurate  for  demonstrating  the 
simulation  procedure. 

For  convenience,  the  differential  equations  of  the  thermal  models  are 
obtained  from  thermal-to-electrical  analogs.  Thermal  capacitance  of 
each  mass  is  the  product  of  that  mass  and  its  specific  heat.  Thermal 
resistance  is  a measure  of  the  temperature  differential  required  to 
obtain  a given  rate  of  heat  flow.  The  thermal-to-electrical  analog 
equivalents  are  given  in  Table  13. 

Table  13 

Thermal-to-Electric  Analogs 


Temperature 

T 

Voltage 

V 

Thermal  Resistance 

1/  Ah 

Resistance 

R 

Thermal  Mass 

mC 

P 

Capacitance 

C 

Heat  Rate 

Q 

Current 

I 

Thermal  Energy 

Q 

Charge 

Q 
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The  electric  circuit  that  represents  the  analog  equations  of  the  thermal 
model  for  the  dc  traction  motor  is  shown  in  Fig.  15.  The  current 
sources  are  the  heat  rates,  i.  e.  , the  power  losses,  of  the  rotor  and 
the  stator.  The  temperatures  of  the  rotor,  air  gap,  and  stator  are 
the  node  variables  with  respect  to  ground,  i.  e.  , the  ambient  tempera- 
ture. The  same  thermal  model  of  Fig.  15,  with  suitable  symbols, 
represents  the  equations  for  the  thermal  model  of  the  alternator. 

The  parameters  required  to  construct  a thermal  model  of  a rotating 
electric  machine  are  not  provided  by  the  manufacturer.  They  must 
be  estimated  from  the  weight  of  the  machine,  the  losses,  and  the 
temperature  when  the  machine  is  operated  at  rated  conditions.  The 
weight  must  be  divided  between  the  rotor  and  stator;  an  average 
specific  heat,  such  as  0.  11  cal/  g/  °C,  is  used  for  the  copper  and  iron 
portions  to  calculate  the  thermal  masses.  The  losses  must  also  be 
divided  between  the  rotor  and  stator.  Thermal  resistances  are  calcu- 
lated by  using  the  rated  temperature,  the  ambient  temperature,  and  an 
estimate  of  the  intermediate  air-gap  temperature.  The  calculations 
for  the  thermal  models  of  the  baseline  electric  machines  are  given  in 


Appendix  A. 
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Thermal  Circuit  Model  for  DC  Traction  Motor 
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Although  the  power  loss  models  provide  five  loss  components  for  the 
alternator  and  the  dc  traction  motor,  the  thermal  models  require  as 
inputs  only  the  total  rotor  loss  and  the  total  stator  loss.  These  losses 
are  allocated  in  the  thermal  models  by  the  two  following  equations, 
which  incorporate  the  weighting  coefficients  ^ and  : 

rdc  dl  1 dZ  Z d3  3 d4  4 d5  5 


P , - S D + S D + S D + S D + S D 

sdc  dl  1 dZ  Z d3  3 d4  4 d5  5 


where: 


R + S = 1.  0 for  i = 1,  2,  ...  5. 
di  di 


In  the  alternator,  all  of  the  armature  power  loss  D occurs  in  the 

stator:  hence,  R = 0 and  S , = 1.  0.  In  the  dc  traction  motor,  the 

dl  dl 

armature  power  loss  occurs  in  the  series  field  and  interpoles  on  the 

stator  and  in  the  armature  winding  on  the  rotor;  the  coefficients  are 

thus  allocated  R ,,  = 0.  95  and  S ,,  = 0.  05.  The  allocations  either  are 
dl  dl 

obvious  or  must  be  made  on  a judgmental  basis.  The  allocations  for 
the  alternator  and  dc  traction  motor  are  given  in  Table  14.  No  matter 
how  the  loss  components  are  allocated,  the  total  power  loss  must  be 
discharged  into  the  total  thermal  mass  of  each  electric  machine. 
Temperature  errors,  if  any,  between  the  rotor  and  stator  of  the  two-mass 
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thermal  model  will  occur  only  because  the  power  loss  components 
are  not  accurately  assigned  to  the  rotor  or  the  stator. 


Table  14 

Weighting  Coefficients  for  Thermal  Models 


of  Electric  Machines 




Mach.i  DC  Machine 

AC  Machine 

me-  n 

ment  r^i 

°2 

°3 

°4 

D5 

Ai 

A2 

A3 

A4 

A5 

Rotor  jjo.  95 

0 

0.6 

0.8 

1.0 

0 

1.0 

0.2 

0.7 

0.2 

Stator|j0.05 

1.0 

0.4 

0.2 

0 

1.0 

0 

0.8 

0.3 

0.8 

The  flow  chart  for  calculating  the  temperatures  of  the  dc  traction 

motor  is  shown  in  Fig.  16.  The  flow  chart  uses  the  parameters  of 

thermal  resistance  R and  R , as  well  as  the  thermal  time  constants 

r s 

R C and  R C . The  rotor  and  stator  components  of  the  power  loss 
r r s s 

are  the  inputs;  the  dc  traction  motor  temperatures  are  the  outputs. 

In  the  dc  traction  motor,  the  heat  from  the  power  losses  is  dissipated 
by  convection  to  the  air,  which  is  blown  through  the  air  gap  by  fans 

3 

at  400  ft  / min.  In  the  alternator,  the  heat  is  dissipated  by  forcing 
an  oil  spray  to  circulate  within  the  machine  at  a rate  of  4 gal / min. 

The  oil  is  cooled  in  an  external  heat  exchanger.  The  numerical  values 


d, 

03 

60 
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Flowchart  of  Thermal  Computations  for  DC  Traction  Motor 
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of  the  parameters  for  the  thermal  models  of  the  electric  machines 
are  given  in  Table  15. 

Table  15 


Thermal  Resistances  and  Time  Constants 
for  Electric  Machines 


DC  Traction 

Parameter  Motor  Alternator 


Rotor-to-gap  therm,  res.  R^ 
Stator-to-gap  therm,  res. 
Gap-to-amb.  therm,  res.  R 

g 

Rotor  therm,  time  constant  t 

r 

Stator  therm,  time  constant  t 

s 


8.  2°  C/kW 
21.  6°  C/kW 
3.  0°  C/kW 
512  s 
1127  s 


18.  8°  C/kW 
7.  4°  C/kW 
3.  0°  C/kW 
65  s 
40  s 


Solid-State  Controller.  The  six  thyristors  of  the  controller  are 
mounted  on  a common  heat  sink  assembly  which  is  cooled  by  convec- 
tion and/  or  forced  air.  The  thermal  model  for  this  component  uses  two 
lumped  thermal  masses:  one  for  the  six  thyristor  junctions;  the  second 
for  heat  sink  assembly.  The  mass  of  the  junctions  is  negligible  compared 
to  the  heat  sink  assembly  and  can  be  neglected.  The  circuit  in  Fig.  17 
represents  the  analog  equations  of  the  controller's  thermal  model. 


thyristor  junction 


Fig.  17 


Thermal  Circuit  Model  for  Solid-State  Controller 
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The  power  loss  is  calculated  as  the  sum  of  junction  losses  of  all 
six  thyristors.  The  junction-to-heat  sink  thermal  resistance  repre- 
sents six  thyristors  in  parallel.  Manufacturers  provide  values  of 
thermal  resistance  for  their  thyristors,  and  data  from  which 
thermal  resistances  and  time  constants  of  their  heat  sinks  can  be 
calculated.  The  parameters  for  the  solid-state  controller  are  given 
in  Table  1 6. 

Table  1 6 


Thermal  Resistances  and  Time  Constants 


for  Solid-State  Controller 


Junction-to-heat  sink  resistance 

+-> 

CD 

58.  8°  C/kW 

Heat 

sink-to-ambient  resistance 

R6 

19.  6°  C/kW 

Heat 

sink  time  constant 

T6 

18.  4 s 
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Changing  Ratings  of  Major  Components 

The  simulation  of  the  electric  propulsion  system  is  intended  for  both 
analysis  and  design  purposes.  The  numerical  values  in  the  models 
of  .the  simulation  are  based  on  the  specific  components  of  the  baseline 
system.  The  simulation  provides  information  on  the  electrical  and 
mechanical  variables,  the  power  losses  and  the  temperatures,  as 
the  vehicle  carries  out  its  mission.  The  specific  components  of  the 
baseline  system  may  have  to  be  changed  several  times  in  the  models 
in  order  to  optimize  the  electric  propulsion  system.  The  changes  in 
the  components  can  result  in  such  effects  as  reduction  of  energy  losses 
over  a mission  cycle;  change  of  recharging  time;  decrease  of  propul- 
sion system  weight;  and  reduction  of  component  temperatures. 

Changes  can  be  made  in  the  electric  machines  of  speed  and  torque 
ratings,  method  of  cooling  and  maximum  temperature.  Changes  can 
be  made  in  the  controller  by  selecting  different  thyristors  and  heat  sink 
assembly.  The  numerical  values  in  the  electric -machine  models  can 
be  modified  either  by  incorporating  the  values  for  the  new  machines, 
or  by  using  dimensional-analysis  on  the  baseline  values  to  find  the  new 
values.  Since  machine  electrical  and  thermal  parameters  vary  slowly 
with  changes  in  power  rating,  the  dimensional-analysis  technique  to  be 
described  in  this  section  is  sufficiently  accurate  for  changes  of  ± 50% 
in  power  rating  from  the  baseline  conditions.  For  small  changes  of 
power  rating,  the  per-unit  machine  ratings  can  be  left  unchanged. 
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Changes  in  Electrical  Ratings  of  Electric  Machine.  The  frame  of  an 
electric  machine  determines  its  torque  rating.  Its  power  rating  is 
proportional  to  how  fast  the  machine  is  run.  Of  course,  many  secon- 
dary factors  must  be  considered,  e.  g.  , mechanical  strength,  losses, 
ventilation,  commutation,  and  bearings.  In  this  section,  we  will 
consider  changes  in  the  power  ratings  of  the  dc  traction  motor  and 
the  alternator,  assuming  that  their  rated  speed  is  not  changed. 

The  power  rating  P of  any  rotating  electric  machine,  such  as  the 
alternator  or  dc  traction  motor,  is  of  the  dimensional  form: 

P = f4JB  N, 
m 

where: 

i = linpar  dimension 
J = current  density 
B^=  flux  density 
N = speed  t 

Assume  that  the  machine  size  is  increased  or  decreased,  i.  e.  , by  its 
linear  dimension  £ t without  changing  either  the  configuration,  the  densi- 
ties, or  the  turns.  Then,  the  copper  and  iron  losses  are  proportional 
3 

to  l , which  is  the  volume  of  material.  The  winding  resistances  both 
in  ohms  and  in  per-unit  terms  vary  as  ^ the  winding  inductances 
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in  henries  vary  as  l , and  the  time  constants  vary  as  a . The  pu 

alternator  reactances  at  a given  speed  and  frequency  vary  as  S.  . The 

friction  and  windage  loss  considered  as  a surface  phenomenon  can 

2 

be  assumed  to  vary  as  I 

Consider  the  70-hp  dc  traction  motor  of  the  baseline  system  as  an 
example  of  how  the  machine  parameters  vary  with  the  linear  dimen- 
sion H . Assume  that  the  dc  traction  motor  is  found  to  operate  with 
too  low  a temperature  over  a mission  cycle.  To  increase  the  temper- 
ature the  rating  will  be  reduced  to  60  hp  and  the  mission  cycle  repeated. 
The  parameters  will  change  as  follows: 


Rating  reduction 

P 

= (60/  70)  = 0.  857 

1/4 

Linear  dimension 

£ 

= (60/  70)  = 0.  962 

Copper  and  iron  loss  (kW) 

P 

= (0.  962)3  = 0.  890 

Resistance  (pu) 

r1 

= (0.  962)  ” 1 - 1.  04 

2 

Friction  & windage  loss  (kW)^ 

= (0.  962)2  - 0.  927  . 

These  parameter  changes  are  small.  For  a 0.  14  pu  reduction  of  power 
rating,  the  pu  resistance  increases  by  only  0.  04  pu;  the  pu  copper 
and  iron  losses  at  the  new  rating  increase  by  only  (0.  890/  0.  85  7)  - 1 


= 0.  04  pu,  as  well. 
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The  assumption  of  maintaining  the  same  number  of  turns  in  the 
electric  machine  as  its  size  is  changed  is  merely  a convenience. 
The  assumption  washes  out  when  the  effects  are  expressed  in  pu. 
After  its  size  has  been  changed  the  machine  can  be  wound  for  any 
combination  of  voltage  and  current  that  meets  the  new  power  rating. 


Changes  in  Speed  Ratings  of  Electric  Machines.  When  the  speed, 
current  or  voltage  ratings  of  an  electric  machine,  such  as  the  alter- 
nator or  dc  traction  motor,  are  changed,  with  or  without  a simul- 
taneous change  in  linear  dimension,  the  change  in  the  parameters  is 
more  complicated  than  for  a change  of  i-  alone.  Under  the  assumption 
of  the  same  winding  geometry  (i.  e.  number  of  slots,  space  factor,  and 
winding  design)  the  voltage  and  current  ratings  can  both  be  varied  by 
a ratio  k,  which  does  not  affect  the  power  rating.  The  ratio  k is 
defined  as  the  turns  ratio  or  the  number  of  turns  of  the  new  machine 
divided  by  the  number  of  turns  on  the  original  machine.  Independent 
values  of  the  turns  ratio  k can  be  assigned  to  the  stator  windings  and 
the  rotor  windings.  The  parameters  are  now  scaled  as  follows: 


Power  rating 
Current  rating 


*4n 

.eV1 


i 2Nk 


Voltage  rating 


1 2 

Resistance  (ft)  ^ k 

2 

Inductance  (H)  £ k 


Reactance  (ft) 

i N k2 

Base  impedance  (ft) 

N k2 

Resistance  (pu)  (£  N) 


Reactance  (pu) 

£ 

Copper  loss  (kW) 

i3 

Copper  loss  (pu) 

(IN)'1 

Iron  loss  (kW) 

iV 

Iron  loss  (pu) 

N £ " 1 

Friction  & windage  loss  (kW) 

£ 2 N3 

Friction  & windage  loss  (pu) 

N2£  "2 

Excitation  power  loss  (kW) 

3 

£ 

Excitation  power  loss  (pu) 

(£  N )~\ 

Consider  the  flywheel  and  alternator  of  the  baseline  system  as  an 
example  of  how  the  machine  parameters  vary  with  the  linear  dimen- 
sion and  the  speed.  Assume  that  the  flywheel  and  alternator  will  be 
changed:  the  rated  speed  will  be  reduced  from  12,000  r/ min  to 
9,000  r/min  and  the  power  rating  raised  from  75  kVA  to  100  kVA. 


The  key  parameters  will  change  as  follows: 
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Rating  increase 

II 

£ 

(100/  75) 

= 1.  33 

Speed  ratio 

N 

(9000/  12,  000) 

= 0.  75 

Linear  dimension 

i 

1 / 4 

(1  . 33/  0.  75)  ' 

= 1.  155 

Resistance  (pu) 

(iN)"1  = 

(1.  155x0.  75)~ 1 

= 1.  155 

Reactance  (pu) 

i 

(1.  155) 

= 1. 155 

Copper  loss  (pu) 

(fN)*1  = 

(1.  55x0.  75)  _1 

= 1. 155 

Iron  loss  (pu) 

Nl  _1  = 

0.  75  x 1.155_1 

= 0.  65 

Friction  & windage 

loss  (pu) 

2 -2 
N i 

0.75^x1.155 

= 0.  42  . 

Changes  in  Electrical  Rating  of  Solid-State  Controller.  The  snubber 
network  loss  varies  directly  with  the  rating  of  the  dc  traction  motor. 

The  thyristor  forward  drop  power  loss  varies  directly  as  the  rated  or 
base  armature  current  of  the  dc  traction  motor,  i.  e.  , the  loss  coefficient 
is: 


R. 


1000 


I 


% 


kW  . 


Changes  in  Thermal  Parameters.  The  parameters  for  the  thermal 
models  of  the  machines  must  be  changed  as  the  ratings  are  altered 
either  by  size  or  by  speed,  or  both.  Two  parameters  are  used  in 
the  thermal  models:  thermal  resistance  and  thermal  mass.  The 
thermal  mass  can  be  combined  with  the  thermal  resistance  to  form 
an  alternative  parameter:  thermal  time  constant.  The  thermal 
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resistances  are  defined  for  the  paths  extending  from  the  rotor  and 
stator  hot  spots  into  the  air  gap,  and  from  the  air  gap  to  ambient - 
temperature  heat  sink.  The  resistances  will  depend  somewhat  upon 
the  speed  of  the  rotor  as  it  affects  heat  transfer  in  the  air  gap.  For 
the  first-order  model,  we  will  assume  that  the  thermal  resistances 
are  independent  of  the  speed.  The  parameters  are  then  as  follows: 

.3 

Thermal  mass  ~£ 

-1 

Thermal  resistance  1 

> 

Thermal  time  constant  r . 


For  example,  for  the  change  of  alternator  rating  from  75  kVA  to 

100  kVA,  and  speed  from  12,000  r/  min  to  9000  r/  min,  in  the  previous 

section,  the  thermal  parameters  will  change  as  follows: 


Linear  dimension 
Speed  ratio 
Thermal  mass 
Thermal  resistance 
Thermal  time  constant 


£ 

N 

3 3 

£ = (1.  155) 

£ "I 

£ 2 =(1.1  55) 2 


= 1 . 1 55 
= 0.  75 
= 1 . 542 
= 0.  867 
= 1. 335 


For  the  solid-state  controller,  as  the  rating  of  the  dc  traction  motor 


current  (I  ) varies,  the  thermal  parameters  are  varied  as  follows: 

% 

Thermal  mass  (C-)  ~ I 

V, 

Thermal  resistance  (R,  and  R,  ) ~ I 

6 6t  % 

Thermal  time  constant  (t  ) — no  change. 
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5.  RESULTS 


Sample  Run 

Using  the  baseline  components,  the  computer  program  for  the  electric 
propulsion  system  simulation  was  run  in  order  to  check  on  several 
items,  namely,  the  models,  the  ratings  selected  for  the  components, 
and  to  validate  the  simulation.  The  results  were  obtained  as  a print- 
out and  as  charts;  the  charts  displayed  traces  of  power,  electrical 
variables,  and  temperatures,  with  time  as  the  independent  variable. 
These  charts  are  shown  in  Figs.  18  through  21. 

Mission  Profile  - The  route  traversed  by  the  existing  Morgantown 
vehicle  was  used  in  the  simulation  to  test  the  baseline  flywheel  energy- 
storage  vehicle.  The  required  elevation  and  velocity  profiles  were 
shown  in  Fig.  4;  they  are  repeated  in  terms  of  roadway  grade  and  ve- 
hicle speed  profiles  in  Fig.  18.  The  route  has  an  outbound  portion, 
which  is  basically  uphill,  to  the  Engineering  Station  (B),  and  a return 
inbound  portion  down  the  hill  to  the  starting  point  at  Walnut  Station  (A). 
The  dwell-times  between  the  segments  of  the  route  are  used  for  charg- 
ing the  flywheel.  The  vehicle  is  assumed  to  make  continuous  round- 
trips,  stopping  at  the  terminals  only  to  handle  passengers  and  charge 
the  flywheel.  The  two  profiles  represent  inputs  to  the  program.  The 
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Vehicle  Speed,  Roadway  Grade  Profile  and  Vehicle 
Tractive  Force  Requirements 


VEHICLE 

SPEED 

(MPH) 
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vehicle  speed  is  shown  as  a reference  on  all  of  the  charts  describing 
vehicle  performance. 

Tractive  Force  - For  vehicle  weight,  wind  conditions,  and  tire 
inflation,  the  computer  calculates  the  tractive  force  required  to  meet 
the  speed  and  grade  profiles.  The  tractive  force  is  shown  in  Fig.  18 
for  the  two  segments  of  the  route.  Negative  tractive  force  obviously 
means  braking.  The  base  tractive  force  for  the  70 -hp  dc  traction 
motor  is  750  lb.  At  620  s on  the  chart,  the  peak  tractive  force  is 
1900  lb;  the  motor  must  deliver  2.  5 times  its  base  value. 

The  peak  tractive  force  would  be  reduced  in  an  actual  vehicle, 
having  an  on-board  speed-control  unit.  The  unit  would  be  designed 
to  control  the  firing  angle  of  the  solid-state  controller  to  coordinate 
as  closely  as  possible  the  actual  speed  of  the  vehicle  with  the  command 
speed.  Furthermore,  the  controller  would  also  limit  the  maximum 
current  that  it  could  deliver  to  the  armature  of  the  traction  motor. 

This  limiter  would  protect  the  controller  and  the  dc  traction  motor, 
as  well  as  limiting  vehicle  acceleration  for  the  comfort  of  the  passen- 
gers. As  a consequence  of  the  control  actions,  the  dc  traction  motor 
would  apply  tractive -force  peaks  less  than  those  shown  in  Fig.  18. 
Instead,  the  vehicle  would  slow  down  while  negotiating  the  maximum 
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grades  and  require  less  tractive  force.  With  reduced  peak  tractive 
forces  and  currents,  the  actual  propulsion-system  losses  would  be 
less  than  the  losses  calculated  in  this  simulation. 

In  the  second  chart  of  Fig.  19,  the  trace  is  the  flywheel  power. 

Negative  values  represent  the  charging  power  when  the  alternator 
operates  a synchronous  motor;  positive  values  represent  the  flywheel 
power  delivered  to  the  alternator  for  propulsion.  In  the  third  chart, 
the  negative  values  represent  the  wayside  inverter  charging  power 
during  the  charging  periods.  During  the  traveling  periods,  the  posi- 
tive values  represent  propulsion  power  and  the  negative  values  repre- 
sent braking  power  at  the  wheels.  The  charging  power  reaches  200  kW 
when  the  flywheel  reaches  full  speed  of  12,000  r/  min.  The  only  regen- 
eration of  energy  to  the  flywheel  in  Fig.  19  occurs  at  about  1000  s,  as 
the  vehicle  is  descending  a 10  percent  grade.  The  braking  power  at 
the  wheels  is  43  kW;  only  7 kW  reaches  the  flywheel.  Of  the  remaining 
36  kW,  6 kW  is  auxiliary  power  and  the  rest  is  the  loss  of  the  electric 
propulsion  system. 

The  power  losses  are  shown  in  the  fourth  chart  of  Fig.  19.  During 
the  charging  period  the  trace  is  the  sum  of  the  power  loss  of  the 
alternator  operating  as  a synchronous  motor  and  the  auxiliary  power. 
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During  the  traveling  period,  the  trace  is  the  total  power  loss  of  the 
electric  propulsion  system  and  the  auxiliary  power.  The  most 
severe  losses  are  evident  when  peak  tractive  force  is  required  on 
the  outbound  segment  of  the  route.  For  example,  at  620  s,  the 
flywheel  power  is  195  kW;  the  propulsion  power  is  85  kW;  and  the 
sum  of  the  losses  and  auxiliary  power  is  1 1 0 kW.  At  this  point  in 
the  mission,  the  propulsion  system  efficiency  is  only  45  percent. 

As  a matter  of  fact,  the  power  losses  are  so  high  on  the  outbound 
segment  that  the  energy  efficiency  of  the  electric  propulsion  system 
is  only  47.4  percent. 

The  trace  of  flywheel  energy  is  shown  in  Fig.  20.  During  the  first 
charging  period,  the  flywheel  is  accelerated  from  half  speed,  i.  e.  , 
25  percent  energy,  to  full  speed,  i.  e.  , 100  percent  energy.  At  the 
end  of  the  uphill  outbound  segment,  the  flywheel  speed  has  declined 
to  7,  760  r/  min,  i.  e.  , 42  percent  energy.  At  the  end  of  the  downhill 
inbound  segment,  the  flywheel  speed  has  declined  to  10,400  r / min, 
i.  e.  , 69  percent  energy.  Note  that  the  7 kW  of  power  regeneration 
at  about  1000  s produces  the  only  increase  of  the  flywheel  energy, 
during  either  the  outbound  or  inbound  segments  of  the  mission  cycle. 

Electrical  Operation  - In  Fig.  20,  the  firing  angle  of  the  solid- 


state  controller  and  the  armature  current  of  the  dc  traction  motor 
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show  the  state  of  the  electrical  system  at  each  interval  of  travel  time. 
Consider  as  a reference  that,  at  a firing  angle  of  90°,  the  controller 
applies  zero  voltage  to  the  dc  traction  motor.  As  the  angle  is  shifted 
toward  0°,  the  voltage  increases  and  the  motor  provides  tractive 
force;  as  the  angle  shifts  toward  180°,  the  voltage  increases  nega- 
tively and  the  motor  provides  regenerative  braking  force.  The  regen- 
erative action  can  be  seen  during  the  downhill  inbound  segment  of  the 
mission  when  the  firng  angle  shifts  to  120°.  At  each  interval  of 
travel  time,  the  computer  reads  out  the  firng  angle  which  makes  the 
propulsion  system  meet  the  speed  and  tractive -force  requirements 
of  the  vehicle.  As  previously  noted,  in  an  actual  vehicle  the  control 
unit  would  generate  the  necessary  firing  angles  as  the  vehicle  traversed 
its  route. 

The  dc  armature  current  trace  in  Fig.  20  roughly  tracks  the  tractive 
force  shown  in  Fig.  18.  Saturation  of  the  dc  traction  motor  causes 
the  armature  current  to  exceed  proportionality  to  tractive  force,  while 
the  compounding  effect  of  the  series  field  acts  on  the  armature  current 
in  the  opposite  direction.  The  peak  armature  current  needed  to  deliver 
the  peak  tractive  force  is  about  2.2  pu  at  about  620  s.  At  this  current, 
the  armature  power  loss  is  about  five  times  its  value  at  rated  current. 
However,  since  the  armature  current  is  less  than  rated  value  over  most 
of  the  route,  the  average  armature  power  loss  is  less  than  the  rated  value. 
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Component  Temperature  - The  temperature  rise  traces  of  the 
critical  parts  of  the  dc  traction  motor,  the  alternator,  and  the  solid- 
state  controller  are  shown  in  Fig.  21.  The  armature  of  the  traction 
motor  reaches  a temperature  of  32°  C near  the  end  of  the  first  uphill 
outbound  segment  and  47°  C near  the  end  of  the  second  uphill  outbound 
segment.  The  rated  temperature  rise  of  the  dc  traction  motor  is  80°  C 
over  a 40°  C ambient.  Since  the  traction  motor  armature  has  a long 
thermal  time  constant  of  512  s,  the  temperature  does  not  track  the 
armature  current.  Furthermore,  the  dc  traction  motor  tends  to  cool 
down  during  the  charging  periods,  when  the  vehicle  is  stationary,  with 
the  result  that  the  dc  traction  motor  temperature  does  not  reach  the 
limits  set  by  the  manufacturer. 

The  stator  of  the  alternator  has  a relatively  short  thermal  time  constant 
of  40  s.  The  alternator  heats  up  during  both  the  charging  and  the  trav- 
eling modes  and  carries  overload  currents  as  well.  The  chart  shows 
a stator  temperature  rise  of  246° C at  the  end  of  the  first  charging 
period,  370° C at  the  end  of  the  first  uphill  outbound  segment,  and 
355° C at  the  end  of  the  second  uphill  outbound  segment.  The  rated 
temperature  rise  set  by  the  manufacturer  is  160°C  over  a 40° C ambient. 
The  alternator  clearly  must  be  replaced  with  a larger  machine,  so  as 
to  reduce  the  power  losses  and  the  temperature  rise. 
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Because  the  junction  thermal  mass  and  time  constant  are  negligible, 
the  trace  of  thyristor  junction  temperature  follows  the  armature  current. 
The  heat  sink  has  a thermal  time  constant  of  only  19  s.  Its  tempera- 
ture rise  reaches  a peak  of  155°C  when  peak  tractive  force  is  required. 
The  rated  temperature  of  the  junction  set  by  the  manufacturer  is  either 
125°C  or  150°C,  depending  upon  the  type  of  thyristor.  For  a 40° C 
ambient,  the  corresponding  rated  temperature  rise  is  either  85° C or 
110°C,  compared  to  the  peak  value  on  the  chart  of  155°C.  The  results 
show  that  the  thyristors  must  be  replaced  with  larger  units  and/  or 
the  thermal  time  constant  of  the  heat  sink  increased  to  insure  that  the 
controller  will  not  fail. 

System  Redesign.  The  purpose  of  testing  the  baseline  electric  propul- 
sion system  was  to  demonstrate  the  procedures  for  using  the  simula- 
tion. The  particular  type  of  propulsion  system  modeled  has  been 
proposed  by  several  companies  for  application  to  a flywheel  energy- 
storage  vehicle.  The  solid-state  controller  and  dc  traction  motor  are 
actually  on  the  existing  Morgantown  vehicle.  The  75-kVA  alternator 
and  the  flywheel  correspond  to  the  energy-storage  unit  originally 
selected  for  the  San  Francisco  MUNI  flywheel  trolley  coach.  The  test 
run  for  the  simulation  of  our  particular  configuration  and  components 
shows  that  the  baseline  dc  traction  motor  is  conservatively  loaded. 
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but  the  baseline  alternator  and  solid-state  controller  are  overloaded. 
In  actual  operation,  they  will  suffer  excessive  temperature  rises  and 
may  be  expected  to  break  down  as  a consequence. 

If  the  baseline  system  were  to  be  built,  one  or  more  of  the  following 
steps  of  redesign  would  be  required  to  insure  that  the  component 
temperatures  did  not  exceed  the  manufacturers'  limits: 

Select  an  alternator  and  a controller  for  larger  power 
ratings  in  order  to  reduce  their  temperature  rises  to 
allowable  levels. 

Reduce  the  size  and  hp  rating  of  the  dc  traction  motor. 
Set  the  base  speed  lower  to  correspond  to  a vehicle 
speed  of  23  mi/h,  at  which  peak  torque  is  required. 

Use  independent  field  weakening  to  reach  the  maximum 
vehicle  speed  of  30  mi/h,  at  which  less  than  peak  torque 
is  required. 

Supply  the  on-board  auxiliary  load  from  a battery,  which 
can  be  charged  from  wayside  equipment  each  time  the 
flywheel  is  charged.  Maintain  the  charge  on  the  battery 
from  the  alternator;  only  allow  it  to  charge  when  it  is 
not  heavily  loaded  with  propulsion  power. 
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Control  the  dc  traction  motor  voltage  with  a combination 
of  alternator  field  control  and  controller  firing  angle 
control  so  as  to  reflect  as  high  a power  factor  as  pos- 
sible to  the  alternator.  The  reduced  armature  and  field 
currents  will  reduce  the  power  losses  and  temperatures 
in  the  alternator'  field  and  stator  windings. 

Increase  the  flywheel  energy  storage  capacity  in  order 
to  reduce  the  speed  drop  of  the  alternator  over  a mis- 
sion. The  alternator  field  current  to  support  the  ter- 
minal voltage  at  large  armature  current  and  low  speed 
will  be  reduced.  The  temperature  rise,  consequently, 

will  be  reduced,  as  well. 

Increase  the  alternator  terminal  voltage  during  the 
charging  period,  so  that  the  voltage  at  the  end  of  the 
period  is  greater  than  1.  2 pu.  For  the  same  current, 
the  power  will  be  greater,  the  charging  period  shorter, 
and  the  stator  temperature  rise  reduced. 

The  proposed  steps  of  redesign  can  be  tested  with  the  simulation  of  the 
electric  propulsion  system  by  making  appropriate  changes  in  the  com- 
ponent models.  The  charts  of  component  temperatures  will  show  when 


the  redesign  has  achieved  a viable  system. 
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Extension  to  Other  Systems 

The  procedure  in  this  report  for  developing  a simulation  of  an  electric 
propulsion  system  for  a flywheel  energy-storage  vehicle  can  be  ap- 
plied to  other  propulsion  systems.  The  simulation  can  be  used  to  ana- 
lyze these  other  propulsion  systems  when  the  components  are  known, 
or  to  select  components  to  meet  criteria  of  temperature  or  power 
losses.  Simulations  are  particularly  suitable  for  determining  the  ca- 
pacity of  the  energy-storage  equipment  required  for  the  vehicle  to 
complete  its  mission. 


In  order  to  develop  a simulation  for  a new  electric  propulsion  system, 
the  models  must  be  formulated  in  the  same  sequence  as  described  in 
this  report.  The  dynamic  models,  which  are  basically  the  equations 
for  the  electric  machines  and  other  components,  will  provide  all  of  the 
dependent  electrical  and  mechanical  variables  of  the  propulsion  system. 
The  power  loss  models  of  the  components  will  then  provide  the  power 
loss  components  as  a function  of  the  electrical  and  mechanical  vari- 
ables. Finally,  the  thermal  models  will  generate  the  component  temp- 
eratures, as  a function  of  the  same  variables. 

For  each  electric  propulsion  system,  the  parameters  for  formulating 
the  models  of  the  components  must  be  obtained  from  a combination  of 
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manufacturers'  data,  handbook  data  of  typical  components,  and  electric  - 
machine  design  texts.  Dimensional  analysis  can  be  used  to  find  param- 
eters for  specific  component  ratings  from  available  data. 

Conclusions 

The  Flywheel  Propulsion  Simulation  has  been  developed  for  use  as  an 
engineering  tool  for  evaluation  and  design  of  flywheel-electric -drive 
propulsion  systems  for  short-range  vehicles.  This  simulation  makes  it 
possible  to  evaluate  performance  of  candidate  propulsion  systems, 
operating  over  given  mission  profiles,  defined  in  terms  of  velocity  and 
grade  profiles.  Wayside  recharge  and  regenerative  braking  are  featured, 
and  power  losses  and  equipment  temperature  rises  are  computed  on  a 
dynamic  basis. 

The  example  shown  in  the  sample  run  (Chapter  V and  Appendix  D)  was 
not  intended  to  be  an  optimum  design,  but  an  illustration  of  how  the 
simulation  can  be  used  in  a specific  instance. 

The  procedure  described  in  this  report  is  a powerful  tool  in  analyzing 
and  designing  electric  propulsion  systems.  The  use  of  this  procedure 
will  help  avert  design  errors,  where  components  are  too  large  or  too 
small,  run  too  hot  or  too  cold,  and  require  expensive  and  time- 
consuming  modifications  after  the  vehicles  are  built  and  tested. 


' 


6.  CONCLUSIONS 


This  report  has  described  a digital  computer  simulation  suitable  for 
modeling  the  electric -drive  system  associated  with  a flywheel  energy- 
storage  vehicle  in  which  the  electrical,  mechanical,  and  thermal  opera- 
tional characteristics  of  the  propulsion  system  are  included.  Parameters 
associated  with  the  route,  the  vehicle,  or  the  propulsion  system  com- 
ponents, can  be  varied  and  the  overall  effect  of  these  variations  on  the 
performance  of  the  vehicle  can  be  rapidly  determined.  The  particular 
simulation  describes  the  system  more  accurately  and  provides  more 
meaningful  data  than  other  previously  used  simulations,  where  efficien- 
cies are  assigned  to  each  component  in  the  drive  train  either  from 
estimated  average  values,  or  values  computed  from  a multivariable 
equation.  Since  the  simulation  includes  thermal  models  of  the  electric 
machines  and  solid-state  controller,  component  temperatures  can  be 
tracked,  and  also  components  of  different  ratings  and  designs  can  be 
evaluated.  Furthermore,  the  simulation  method  can  be  applied  to 
electric -drive  systems  for  battery  and  flywheel  energy-storage  vehicles, 
to  hybrid  vehicles,  and  to  vehicles  operating  from  wayside  power. 


For  the  sample  run  and  prototype  drive  system  described  in  this  report, 
the  results  indicate  that  the  alternator  must  be  larger  (higher  rating)  to 
reduce  its  operating  temperatures  to  acceptable  levels.  Likewise,  the 
dc  traction  motor  can  be  smaller,  because  its  operating  temperature  is 
less  than  the  allowable  level.  The  rating  of  the  solid-state  controller 
should  be  increased  to  reduce  its  operating  temperature. 

The  results  of  our  preliminary  work  have  shown  that  the  electric -drive 
system  can  be  optimized  for  low  losses  and  maximum  usage  of  the 
flywheel  stored  energy.  The  following  steps  are  required.  First,  the 
machines  should  be  operated  for  as  little  time  as  possible  beyond  their 
ampere  ratings.  A current  limiter  should  be  incorporated  into  the  con- 
troller, or  larger  machines  should  be  used  if  necessary.  Second,  the 
flywheel  and  alternator  of  the  drive  system  should  not  be  burdened  with 
supplying  the  auxiliary  power  at  all  times , particularly  when  the  drive 
system  is  operating  at  full  capacity.  A battery  should  be  used  as  a 
supplementary  power  source  for  the  auxiliary  load.  Third,  a phase - 
controlled  rectifier  should  not  be  used  to  control  the  speed  of  the  dc 
traction  motor,  particularly  when  the  dc  motor  is  required  to  run  at  low 
speed  and  high  torque  over  a substantial  portion  of  the  route.  The  high 
current  at  low  power  factor  drawn  from  the  alternator  causes  its  stator 
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and  field  losses  to  be  excessive.  To  reduce  system  losses,  the  dc 
motor  speed  should  be  controlled  by  the  alternator  field  current;  the 
controller  should  be  operated  either  full-on,  firing  angle  a = 0°  , or  full 
regeneration,  a'cslSO". 


The  electric -drive  system  described  in  this  report  is  frequently  used 
for  studies  of  flywheel  energy-storage  vehicles.  For  this  reason,  it  was 
selected  as  a basis  for  developing  the  modeling  procedure.  Other  con- 
figurations, which  are  less  lossy  over  a given  route,  can  be  studied  by 
the  same  techniques.  Electric -drive  systems  which  are  patterned  after 
stationary  industrial  drives,  or  after  vehicle  drives  using  contact-rail 
power,  will  not  necessarily  be  optimum  for  flywheel  energy-storage 
vehicles.  The  regenerative  braking  feature  of  the  system  described  in 
this  report  would  probably  be  more  attractive  if  a true  urban  route  (over 
five  stops  per  mile)  were  used  to  test  the  drive  system,  since  braking 
would  be  done  more  frequently. 

In  summary,  this  electric -drive  system  simulation  can  be  used  to 
determine  suitability  of  a given  system  configuration  for  an  intended 
mission.  Also,  changes  in  the  ratings  of  components  and  different 
configurations  of  components  can  be  analyzed  in  order  to  optimize  the 


drive  system. 


- . 
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APPENDIX  A,  MATHEMATICAL  MODELS 
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Part  I 


DYNAMIC  AND  POWER  LOSS  MODELS 


The  dynamic  and  power  loss  mathematical  models  are  the  sets  of 
equations  that  describe  the  flywheel  energy-storage  vehicle's  drive 
system  in  both  its  traveling  and  wayside  charging  modes  of  operation. 
For  convenience,  the  equations  are  divided  into  sections  corresponding 
to  the  components  or  functions  which  they  describe.  A pictorial  diagram 
of  the  major  system  components  is  shown  in  Fig.  3.  The  input  con- 
stants for  the  baseline  drive  system  are  given  in  Appendix  B.  The 
parameters  for  the  equations  are  listed  in  Appendix  C.  The  prime 
objective  of  the  calculations  using  the  dynamic  models  is  to  determine 
the  flywheel  power  at  each  instant  of  time  as  the  vehicle  travels  on  its 
route.  Hence,  the  calculations  start  from  the  tractive  force  and  power 
required  by  the  vehicle  drive  wheels,  and  are  carried  back  through  the 
dynamic  models  of  the  propulsion-system  components  to  the  flywheel. 
The  equations  are  presented  in  this  same  order  of  calculation. 

A.  Propulsion  Requirements 

The  input  variables  from  the  mission  profile  are: 


Vehicle  velocity  V{t)  ft/  s 


Roadway  grade 


G(t)  rad 


Wind  velocity 


mi/  h 


The  tractive  force  components  are: 


Inertia  force 


Grade  force 


F = W G 
g 


Tire  friction  force 


F = n W ( 1 + C V) 
t c 


Aerodynamic  drag  force 

F = i-  p (V+  V )2  SGN  (V  + V ) C A 
w 2 r w w d 


The  total  tractive  force  and  traction  power  at  the  drive -wheels 

Total  tractive  force  F =F.  + F+  F+  F 

l g t w 

-3 

Total  traction  power  = 1.  36  x 10  F*V 


If  F and  are  negative. 


the  dc  traction  motor  must  shift  into 


lb. 


lb. 


lb. 


lb. 


are: 

lb. 


kW 


the 


regenerative  braking  mode. 


A 


B. 


Gearbox  (Transmission) 


The  vehicle  speed  and  tractive  force,  when  the  dc  traction  motor  is 


operating  at  its  base  speed  of  2730  r/  min  and  70.  17  hp,  are  given  by: 

2-ttN  R 

Vehicle  speed  V = — — = 47.  18  ft/  s 


Eq 
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Tractive  force 


5252  7 igb 


N. 


R 


751 


lb 


The  gearbox  losses  for  an  assumed  efficiency  p are: 


Motoring 


kW 


Braking 


1 - 


'gb 


P, 


kW 
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C. DC  Traction  Motor 

The  normalized  vehicle  speed  and  tractive  force  are. 
Vehicle  speed,  normalized  V 

o 


V 

oVb 


pu  1 


Tractive  force,  normalized 

(motoring)  F 

o 

(braking)  F 

o 


2 

pu 
pu  . 


The  equations  for  the  armature  current  needed  to  develop  the  required 
tractive  force  must  take  into  account  the  series  field,  saturation,  and 
field  weakening  control  above  the  base  speed  of  the  dc  traction  motor. 

Ampere  turns,  ratio  of  series-to-total,  at  dc  traction  motor  base 
conditions  for  k se ries -to-shunt  field  turns  and  c field  weakening. 


C 


\ a I 

ab 

I.  + \ a I , 
f , ab 

dcb 


pu 


Armature  current  parameter 


SQ  = 


/ 1 -C 

\ 

- 

, — 

|+  B 

F 1 

j dc 

o| 

pu. 


Armature  current  normalized 


IQ  = SQ  + ^jSQZ  + 


r 


A + B (1  -C  ) 
dc  dc  1 


pu  . 


5 
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When  the  dc  traction  motor  runs  over  its  base  speed  a and  a < 1, 
the  shunt  field  is  assumed  to  be  weakened.  The  normalized  armature 
and  shunt  field  currents  are: 


Armature  current,  normalized 


I 

o 


E 


bb 


V 


dc. 


V 

F 

a ° 


pu 


Eq. 
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1 - C, 


dc 


E 


bb 


V 


V 


-B 


dc. 


dc 


-CI 
] o 


pu  . 


The  normalized  ampere  turns  for  shunt  and  series  fields  are  given  by: 


Full  field  I = ( 1 - C , ) + C , I pu  7 

f , 1 1 o 

dc 


Field  weakened  F = (1-C,)I„+C1I  pu  „ 

f , 1 f 1 o 

dc  r 


The  shunt  and  series  fields  of  the  dc  traction  motor  are  reversed  when 
the  motor  transfers  from  the  motoring  to  the  braking  mode.  The  vol- 
tage equations  are  the  following: 

Discriminant  SGN  F = + 1 FOR  F ^ 0 (motoring)  8 

SGN  F = - 1 FOR  F < 0 (braking)  , 


A -6 


Flux-conductor  product 


<f> 


SGN  F 


dc 


dc 


dc 


B,  I1 , 
dc  f 


dc 


pu 


Eq 
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Back  emf  base  E,  = V ( 1 - ^ - 2.  0 

bb  dc, 
b 


V 


10 


Back  emf  normalized 
E , 


bb 


V 


dc. 


V $ , 

0 dc 


pu  . 


1 1 


Terminal  voltage  normalized 

V 


E 


dc 


o V 


dc. 


E + T I + — 
4 o V 


pu 


dc. 


12 


When  the  base  speed  at  which  field  weakening  is  used  is  reduced 
to  a and  a smaller  traction  motor  is  used,  the  motor  parameters 
change  to: 


Armature  resistance 


R 


R 


dc 


dc 


1 3 


Armature  current,  base 
.1 


I = cr  I 


% 


14 
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The  five  loss  components  and  the  efficiency  for  the  dc  traction  motor 
are  given  in  the  following  equations: 


Armature  copper  loss,  120°C 


D, 

= 1.39  x 10  3 

(I  1 I )2  R 1 
ab  ° adc 

kW  . 

1 5 

Shunt  field  loss 

- ,0'3  V vt 

dc 

kW  , 

1 6 

dc 

Friction,  windage  and 

core  loss 

°3 

2.  5 

= D V 
3b  ° 

kW  . 

17 

Stray  load  loss 

°4 

0.01  -p 

" n , 2 

gb 

(motoring) 

kW 

18 

°4 

= °-01'bgb*|P 

j (braking) 

kW  . 

Brush  loss 

D5 

= 10'3*  2 (I  1 

I ) 
o 

kW  . 

19 

ab 

Total  loss 

P5 

= di  + d2+d 

+ D+D 
3 4 5 

kW  , 

20 

DC  traction  motor  efficiency 

1 

(motoring) 

^dc 

1 + 

P5  11  gb 
P2 

1 

(braking) 

ndc 

1 + 

P5 

"gb  |P2l 

21 


D. 


Solid-State  Controller 


The  controller  losses  and  efficiency  are: 


Snubber  loss  R 


,o'!  (x  inX'%  ■•) 


kW 


Thyristor  loss  R_  = a R~  I 

4 Z,o 

b 


kW 


Total  controller  loss 
P, 


Rl  + R2 


kW  . 


Controller  efficiency 
r> 


dr 


1 + 


100  R 


(motoring)  pu 


dr 


100  R 


(braking)  pu 


to/  from 
alternator 


Fig.  A-l  Circuit  Diagram  for  Solid-State  Controller 


E. Alternator 

The  alternator  power,  reactive  power,  and  kVA  are  given  in  the 
following  sequence  of  equations: 

Power  to  controller 


P = P_  + P . + P_  + P,  - D. 
pdr  Z 4 5 6 < 


kW 


Alternator  power  output 

POW  = P , +P,  kW  . 

pdr  1 


Reactive  power  to  controller,  (fundamental  component) 


rkVa  . 


Auxiliary  load  reactive  power 

Qj  = P tan  0 rkVa  , 


Alternator  reactive  (fundamental  component)  power  output 


Q i 4-  Q 
pdr  1 


A-9 


Eq. 
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2 


3 


4 


Q 


rkVa 


5 


A - 1 0 


Alternator  kVA  output,  (fundamental  component) 
S. 


1 


I 2 ~ ' 

^POW  + Q 


kVA 


:q. 

6 


Alternator  power  factor 


cos  9 = 


POW 


pu 


Alternator  (harmonic  components)  reactive  power  output 

Q = 3 I ll'!  (0.  242)  VI  * I 
h 1 ab° 


rkVa 


Alternator  total  kVA  output, 
S 


V^,! 


kVA 


The  alternator  base  and  normalized  'roltage  and  current  in  terms 
of  the  rewind  factor  RWF  are  given  in  the  following  equations: 


Alternator  base  voltage  ( £ -n) 


V = V (RWF) 
b tb 

o 


Alternator  base  armature  current 


V 


3 

10  SL 
t 

3 V 


10 


1 1 


Alternator  normalized  armature  current 


A- 1 1 


I03s 

4 " 3V,  >4 


Alternator  normalized  fundamental  current 


Eq. 

pu  12 


pu  13 


The  following  equations  describe  the  phasor  diagram  model,  such 
as  shown  in  Fig.  . The  equations  are  used  to  calculate  the  excita- 
tion voltage  E , taking  into  account  saturation. 

5 


Alternator  power  angle 
tan  <5 


I. 

X gj  cos  9 - R sm  0 

4 

q o a 

f 

L ac  J 

V. 


t 


\ + S 


X gj  sin  0 + R cos  9 
q o a 

ac 


pu 
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Alternator  air-gap  voltage 


pu  . 


15 


A- 1 2 


Field  current  for  E 


ac 


A E 
ac  r 

to  - B E 
o ac  r 


pu  . 


Eq. 
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Slope  of  air-gap  line  at  to 


K = to  K 
agr  o ag 


pu 


1 7 


Slope  of  Kingsley  air-gap  line  (cannot  exceed  K 


ag. 


slope  B = 


ac 


pu 
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Ratio  of  slopes  K 


K 


ag. 


fact  slope  B 


pu 


19 


Kingsley  saturated  synchronous  reactance 


X 


= X + 


X , - X 
d a 


sat 


a K 


fact 


pu 


20 


Excitation  voltage 

E - ^ cos  5 + I 
*b  f 


R cos  ( <5  + 0)  + X,  to  sm  (5  + 0) 
a d , o 

ac  sat 

pu 
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Alternator  field  current 
I 


5 slope  B 


pu 


22 


'5  - I5  (pu)  ■ I 


A 


A- 1 3 


The  firing  angle  of  the  controller  thyristors  needed  to  provide  the 
required  dc  traction  motor  armature  voltage  is  given  in  the  following 
equations : 


Commutation  voltage  drop 

1 1 

3 1 

1C  - — — I I X Z oj 

x tt  a.  o d b o 
b 


F-q. 

V . 23 


Controller  firing  angle 

cos  c y 


V , 1C  I 1C  < 3.  0 
dc^  o x 


2.  34  V 

t 


pu  . 24 


The  loss  components  and  efficiency  of  the  alternator  are  given 
in  the  following  equations: 


Armature  copper  loss 


A. 


Exciter  and  field  loss 


A. 


Core  loss 


kW  25 


kW  . 26 


kW  . 27 


A-14 


Friction  and  windage  loss 


A . = A . u> 

4 4,  o 

b 


kW 


Eq. 

28 


Stray  load  loss 


(\ 

AS  A5  + A5  l RWF 

bo  bi  V 


) 


kW  . 29 


Total  alternator  loss 


P7  ‘ Al  + A2+  A3+  A4+  A5 


kW  30 


Alternator  efficiency 


ac 


1 + 


POW 


(traveling)  pu  31 


^ac  1 |POW| 


(charging)  pu 


A- 1 5 


F. Flywheel 

In  the  traveling  mode,  the  flywheel  supplies  the  net  power  for  propulsion, 
losses,  and  auxiliary  load.  The  equations  for  calculating  the  flywheel 
acceleration  and  angular  speed  are  as  follows: 

Total  drive  system  loss  gq 

P„  = P + Pr  + P.  + P„  kW  . 1 

3 4 5 6 7 

Total  flywheel  power 

P = P + P + P - P.  kW  . 2 

12  3m 


Flywheel  integration  constant 


K, 


550 


0.  746  gj 


2 

f 

max. 


1 . 5562  x 1 0_5kW_  1 s~  1 3 


for 


Flywheel  inertia  coefficient 


L = 30  lb-  ft-s 
w = 400-rr  rad/  s 


max. 


L 


K, 


pu 


Flywheel  maximum  speed  coefficient 


pu  . 


5 


A-l  6 


Flywheel  acceleration 


Eq 


'K.  K2  K3  “f  IT 
max  \ o 


rad/  s' 


Flywheel  angular  speed 


CO  = to  + 

o o. 

1 


a 


to„  \ , /to 
f \At  / f 


pu  , 7 


max , 


max; 


A - 1 7 


Part  II THERMAL  MODELS 

The  thermal  models  are  the  equations  that  provide  the  component 
temperature  rises  in  terms  of  the  component  losses  calculated  from 
the  power  loss  models  of  Part  I. 


Eq. 

1 


G. 


DC  Traction  Motor 


Loss  power  to  rotor 


P = R D.+  R,  D + R,  D + R , D + R , Dr 
r,  d.  1 d 2 d.  3 d.  4 d 5 

dc  1 2 3 4 5 


kW 


Loss  power  to  stator 


P = S,  D.+  S,  D + S D + S D + S D 
s , d,  1 d 2 d0  3 d.  4 dc  5 

dc  1 2 3 4 5 


Constraint 


R . + S . = 1.0  for  l = 1 , 2,  ...  5 

d.  d . 

i i 


Rate  of  rotor  temperature  rise 


kW  . 2 


T 


T 


dc 


- T + P R 
2 r r r 

^dc  dc  dc  dc 


dc 


C/s  . 


A-l  8 


Rate  of  stator  temperature  rise 


T 


dc 


- T + P R 
^dc  Sdc  Sdc  Sdc 


dc 


Eq 

'C/s  . 5 


Air-gap  temperature  rise 


T 


R R T + R R T 
g , s,  r,  g,  r,  s. 
&dc  dc  dc  6dc  dc  dc 

t R R+R  R + R R 

dc  r , s , r g . s , g , 

dc  dc  ac  dc  dc  dc 


Loss  power  to  air  gap 


T 


’dc 


dc  R 

out  g 


dc 


Rotor  time  constant 


r = R C 

r r r 

dc  dc  dc 


Stator  time  constant 


r = R C 

s , s , s , 

dc  dc  dc 


kW  , 7 


s . 8 


s . 9 


A-I  9 


H.  Alternator 

Loss  power  to  rotor 

Eq. 


p 

r 

ac 

= R A + R A+R  A+R  A+R  A 
a,  1 a2  2 a3  3 a4  4 &5  5 

kW  . 1 

Loss  power  to  stator 


P 

s 

ac 

= S A + S A + S A + S A + S A 
ai  1 a2  2 a3  3 a4  4 a5  5 

kW  . 2 

Constraint 

R + S =1.0  for  i = 1,  2,  ...  5 3 

a.  a. 

l i 

Rate  of  rotor  temperature  rise 


T 

r 

ac 

T - T + P R 

g r r r 

ac  ac  ac  ac  0~,  . 

T C/s  . 4 

r 

ac 

Rate  of  stator  temperature 


T 

s 

ac 

T - T + P R 

g s s s 

ac  ac  ac  ac  0„,  r 

= C/  s , 5 

s 

ac 

A-20 


Air-gap  temperature  rise 


R R T + R R T 
g s r g r s 

ac  ac  ac  ac  ac  ac 


R R + R R + R R 
r s r g s g 

ac  ac  ac  ac  ac  ac 


Eq. 

6 


Loss  power  to  air  gap 


T 


ac 


out 


R 


ac 


ac 


kW  . 7 


Rotor  time  constant 


r 

ac 


R C 
r r 
ac  ac 


Stator  time  constant 


s . 8 


s 

ac 


R C 
s s 
ac  ac 


s . 9 


A-21 


I.  Solid-State  Controller 


Rate  of  heat  sink  temperature  rise 


Eq. 

°C/s  . 1 


Loss  power  to  ambient  air 


P 


6 

out 


kW  2 


Heat  sink  time  constant 


T 


6 


R C. 
6 6 


s , 3 


Thyristor  junction  temperature  rise 


T 


6 

t 


T6 


t P 


o 


c . 


4 


Part  III  WAYSIDE  CHARGING  MODE 


The  auxiliary  power  (P^)  is  set  to  zero.  The  tractive  power  (P^) 

and  wind  force  (V  ) are  also  zero.  The  vehicle  transmission 
w 

losses  (P  ),  dc  traction  motor  losses  (Pc),  and  controller  losses 
(P^)  are  all  zero.  The  equations  representing  the  inputs  to  the 
alternator  (operating  as  a synchronous  motor)  as  it  accelerates  the 
flywheel  from  initial  speed  to  full  speed  are: 

Terminal  voltage 


V = V * co 

t to 

c 

max 


Power  input 


P. 

in 


w 


= K • P. 
p.  m 

in  max 


CO 

o 


kVA  input 
S 


P. 

in 

w 


| cos  0| 


V . 


kW  . 


kVA  . 


Armature  current 


A-23/A-24 


The  other  alternator  parameters  and  the  alternator  losses  in  the 
charging  mode  are  computed  by  using  the  same  equations  as  those 
used  in  the  traveling  mode,  with  two  exceptions: 

(No  harmonics  in  wayside  power) 

4f 

(p»„  ■ p>) 

r)  = — (Alternator  efficiency  during  charging) 

ac  x . 

in 

w 

where  P is  the  total  alternator  (motor)  loss  in  the  charging  mode. 


Flywheel  acceleration  and  angular  velocity  are  found  from. 


Acceleration 


° - K1  K2  K3 


. P.,  - P. 

/ 7 m \ 


Eq. 

Z 

rad/  s . 5 


max 


Angular  velocity. 


GO 


0 


+ 


’ At 


rad/  s . 6 


APPENDIX  B.  INPUT  CONSTANTS  FOR  COMPUTER 

PROGRAM 

B-  1 

Part 

I CONSTANTS  FOR  A COMPUTER  SIMULATION  RUN 

Symbol 

Symbol 

Nom. 

(Program)  (Equations) 

Definition 

Value 

Units 

WT,  VMASS 

W 

Weight  of  vehicle  with  load 

12,238 

lb 

VAREA 

A 

Projected  area  of  vehicle 

53. 

ft2 

A A 

M- 

Rolling  coefficient  of  friction 

0. 0230 

-- 

B B 

C 

c 

Coulomb  coefficient  of  friction 

0. 000175 

h/  mi 

C C 

cd 

Aerodynamic  drag  coefficient 

0.  610 

-- 

WIND 

V 

w 

Wind  velocity 

30. 

mi/  h 

I 

I 

Flywheel  moment  of  inertia 

45. 

slug-ft 

ALIM 

V 

max 

Acceleration  limit  (abs.  ) 

1 . 0 

mi/  h / s 

F B 

L 

Base  value,  tractive  force 

751 

lb 

V B * 

L 

Base  value,  vehicle  velocity 

47.  1 8 

ft/  s 

IAB 

I 

ab 

Base  value,  dc  motor  armature  current 

140. 

A 

W 4> 

CO 

o 

Speed  of  energy  storage  unit  (initial) 

0.  5 

pu 

NGB 

ngb 

Efficiency  of  gearbox 

0.  92 

-- 

VDCB 

X 

Base  value  of  dc  motor  terminal  voltage 

420. 

V 

VTRUN 

v, 

run 

io5Ki 

Const,  terminal  voltage  ( l-n)  of  alt  .-traveling  215. 

V 

K1 

Integration  constant 

1 . 55626 

-1  -1 

kW  s 

RGDC 

R 

gdc 

DC  motor  thermal  res.  , gap  to  ambient 

2.  97 

°C/  kW 

RADC 

R 

adc 

DC  motor  arm.  circuit  res.  at  20°  C 

0.  1 50 

ohms 

D3B 

°3 

DC  motor,  friction,  windage ,&  core  losses 

at  base  conditions 


B-2 


(continued) 


Symbol 

Symbol 

Nom. 

(Program)  (Equations) 

Definition 

Value 

Units 

AAC 

A 

ac 

AC  machine,  saturation  curve,  Frolich 

coefficient 

0. 3584 

pu 

BAC 

B 

ac 

AC  machine,  saturation  curve,  Frolich 

coefficient 

0.  6417 

pu 

ADC 

Adc 

DC  motor,  saturation  curve,  Frolich 

coefficient 

0.  50 

pu 

BDC 

Bdc 

DC  motor,  saturation  curve,  Frolich 

coefficient 

0.  50 

pu 

PSI 

'I' 

IR  drop/  dc  motor  base  term,  voltage 

0.  05 

pu 

R2B 

\ 

SCR  loss  coefficient,  controller 

0.  412 

kW 

A1B 

\ 

AC  machine,  armature  Cu  loss  at  base  cond 

3.  938 

kW 

A2B</> 

\ 

AC  machine,  exciter  and  field  Cu  loss,  const  0.  28 

kW 

A2B  1 

A ° 
A2 

AC  machine,  and  field  Cu  loss,  coefficient 

0.  2394 

kW 

A3B 

A3b' 

AC  machine,  core  loss,  coefficient 

1 . 601 

kW 

A4B 

A4 

AC  machine,  friction  &.  windage  loss  at 

4b 

maximum  speed 

2.  045 

kW 

A5B  (f> 

A5 

AC  machine,  stray  load  loss,  constant 

0.  207 

kW 

- 

5b 

A5B1 

V 

AC  machine,  stray  load  loss,  coefficient 

1 . 671 

kW 

PI 

P1  1 

Real  power  to  auxiliary  load 

6.  0 

kW 

T81 

9 

1 

Power  factor  angle,  auxiliary  load 

0. 64350 

rad. 

B-3 


(continued) 


Symbol  Symbol 
(Program)  (Equations ) 


XQ 

X 

q 

RAAC 

R 

a 

nac 

X2<£ 

Xd 

XD 

xd 

WFMAX 

Cl) 

f 

max 

LAMBDA 

X 

RRDC 

R 

rdc 

RSDC 

R 

Sdc 

RCRDC 

T 

rdc 

RCSDC 

T 

Sdc 

RD1 

\ 

RD2 

\ 

RD3 

\ 

RD4 

V 

Nom. 


Definition 


Value 


AC  machine,  q-axis  synchronous  reactance  1.24 

AC  machine,  armature  resistance  0.  0495 

AC  machine,  armature  reactance  0.145 

AC  machine,  d-axis  unsaturated  synchronous 
reactance  3.40 

Maximum  speed,  flywheel  & alternator  12,000. 
Series/ shunt  field  turns  0.  04034 

DC  motor,  thermal  resistance,  rotor-to-gap  8.  17 
DC  motor,  thermal  resistance,  stator - 
to-gap  21.6 

DC  motor,  rotor  thermal  time  constant  512. 

DC  motor,  stator  thermal  time  constant  1127. 

Thermal  weighting  coefficient,  -D^,loss 

to  dc  motor  rotor  0.  95 

Thermal  weighting  coefficient,  -D^  loss 

to  dc  motor  rotor  0. 


Thermal  weighting  coefficient,  D^,  loss 
to  dc  motor  rotor 


0.  6 


Thermal  weighting  coefficient,  D^,  loss 


0.  8 


Units 

pu 

pu 

pu 


pu 

r/  min 


°C / kW 


°C / kW 


s 

s 


to  dc  motor  rotor 


B-4 


Symbol  Symbol 
( Program  ) (Equations) 


(continued) 


Definition 


Nom. 

Value  Units 


RD5 

\ 

Thermal  weighting  coefficient,  D , loss 
to  dc  motor  rotor 

1 . 0 

SD1 

% 

Thermal  weighting  coefficient,  , loss 

to  dc  motor  stator 

0.  05 

SD2 

% 

Thermal  weighting  coefficient  - D , loss 
to  dc  motor  stator 

1 . 0 

SD3 

% 

Thermal  weighting  coefficient,  loss 

to  dc  motor  stator 

0.  4 

SD4 

\ 

Thermal  weighting  coefficient,  D^,  loss 
to  dc  motor  stator 

0.  2 

SD5 

% 

Thermal  weighting  coefficient,  loss 

to  dc  motor  stator 

0. 

R6 

R6 

Thermal  resistance  of  controller  heat 
sink  to  ambient 

1 9.  6 

TAU6 

t6 

Controller  heat  sink,  thermal  time  const. 

1 8.  4 

SIGMA 

a 

Field  weakening  factor  new/  old  motor 
base  speed 

1 . 0 

1 FDCB 

Base  value,  shunt  field  current 

12.  0 

TB 

\b 

Base  output  torque,  dc  motor 

1 35.  0 

°C/  kW 


s 


A 

lb"  ft 


B-5 


(continued) 


Symbol 

Symbol 

Nom. 

(Program)  (Equations  ) 

Definition 

Value 

Units 

NB 

N 

b 

Base  speed,  dc  motor 

2730 

rj  min 

HPB 

Hpb 

Base  horsepower,  dc  motor 

70.  1 7 

hp 

VFDC 

vf 

Shunt  field  voltage,  dc  motor 

24.  0 

V 

dc 

RWF 

RWF 

Rewind  factor,  ac  machine,  new/ old 

term,  -volt 

3.  0 

-- 

ER<J 

E 

r 

Air  gap  voltage  at  base  condition 

1. 05 

pu 

RRAC 

o 

R 

AC  machine  thermal  resistance,  rotor- 

rac 

to-gap 

18.  8 

°C/  kW 

RSAC 

R 

AC  machine,  thermal  resistance,  stator- 

ac 

to-gap 

7.  42 

°C / kW 

RGAC 

R 

gac 

AC  machine,  thermal  resistance,  gap-to-amb.3.  03 

°C/  kW 

RCRAC 

T 

AC  machine,  rotor  thermal  time  constant 

64.  9 

s 

rac 

RCSAC 

T 

AC  machine,  stator  thermal  time  constant 

39.  5 

s 

ac 

SB 

sb 

Base  kVA  output  of  alternator 

75.  8 

kVA 

VTB$ 

\ 

AC  machine,  base  term,  voltage  ( l-n) 

o 

before  rewind 

115.0 

V 

KAG 

K 

ag 

Slope  of  air  gap  line 

1.519 

pu 

XA 

X 

a 

AC  machine,  armature  leakage  reactance 

0.  085 

pu 

B-6 


(cont  Lnued) 


Symbol 

Symbol 

Norn. 

(Program  ) 

(Equations  ) 

Definition 

Value 

Units 

I5B 

AC  machine,  field  current,  base  condition 

14.  26 

A 

5b 

RA1 

R 

ai 

Thermal  weighting  coefficient-A  j loss  to 

ac  machine  rotor 

0. 

-- 

R A2 

R 

a2 

Thermal  weighting  coefficient-A^  loss  to 

ac  machine  rotor 

1.0 

RA3 

R 

a„ 

Thermal  weighting  coefficient-A^  loss  to 

3 

ac  machine  rotor 

0.  2 

RA4 

R 

a„ 

Thermal  weighting  coefficient-A^  loss  to 

4 

ac  machine  rotor 

0.  7 

-- 

RA5 

R 

a^ 

Thermal  weighting  coefficient,  A^  loss  to 

5 

ac  machine  rotor 

0.  2 

-- 

SAl 

S 

ai 

Thermal  weighting  coefficient,  A loss  to 

ac  machine  stator 

1. 0 

SA2 

S 

a2 

Thermal  weighting  coefficient,  A^  loss  to 

ac  machine  stator 

0. 

-- 

S A3 

S 

a„ 

Thermal  weighting  coefficient,  A^  loss  to 

3 

ac  machine  stator 

00 

o 

-- 

S A 4 

S 

a4 

Thermal  weighting  coefficient,  A^  loss  to 

ac  stator 


(concluded) 


Symbol 
( Program) 

SA5 

KP  IN 
R6T 

VTCMAX 

PIN  MAX 

MF 

ZB 


Symbol 

(Equations) 


Definition 


Nom. 

Value 


max 


P. 

in 

max 


MF 


Z 


b 


Thermal  weighting  coefficient,  A loss 

D 

to  ac  machine  stator  0.  8 

Input  power  gain  factor,  - charging  mode  1. 0 

Thermal  resistance,  controller  SCR  to 

heat  sink  5 8.  8 

AC  machine, terminal  voltage  (max), 

charging  mode  414.  0 

Maximum  input  power,  ac  machine, 

charging  mode  200.  0 

Mass  factor  due  to  rotary  inertia  1. 08 

Base  impedance  of  alternator  4.  70 


B-8 


Symbol 

Symbol 

Part  II  OTHER  CONSTANTS 

Nom. 

(Program  ) 

(Equations) 

Definition 

Value 

Units 

EBB 

Ebb 

DC  motor,  back  emf  at  base  condition 

397. 

V 

g 

Gravitation  constant 

32.  17 

ft/s2 

V 

Nominal  flywheel  moment  of  inertia 

30. 

slug  -ft2 

R 

Drive  wheel  radius 

1. 1833 

ft 

V 

Gear  ratio  (motor  speed/  wheel  speed) 

7.  17 

-- 

P 

Air  density  (sea  level)  2.  378 

-3 

x 10 

3 

slugs/  ft 

“f 

max 

Nominal  flywheel  maximum  speed 

400  TT 

rad/  s 

D 

00 

o . 

Initial  value  of  flywheel  speed 

0.  5 

pu 

MTRIP 

1 

Number  of  round  trips  for  a run 

1,5 

- 

APPENDIX  C.  PARAMETER  LISTING  FOR  COMPUTER  PROGRAM 


C-I 


Symbol  Symbol 

(Program)  (Equations) 

PARAMETERS  FOR  COMPUTER  PROGRAM 
Definition 

Units 

A 1 

Ai 

AC  machine,  armature  cu  loss,  200°C 

kW 

AZ 

A2 

AC  machine,  exciter  and  field  cu  loss 

kW 

A3 

A3 

AC  machine,  core  loss 

kW 

A4 

A4 

AC  machine,  friction  and  windage  loss 

kW 

A5 

A5 

AC  machine,  stray  load  loss 

kW 

Cl 

C1 

(Series/ total)  amp-turns  at  dc  motor  base 

-- 

D1 

D, 

Armature  circuit  cu  loss  at  120°C 

kW 

DZ 

°2 

DC  motor,  shunt  field  cu  loss 

kW 

D3 

°3 

DC  motor,  friction,  windage,  and  core  losses 

kW 

D4 

°4 

DC  motor,  stray  load  loss 

kW 

D5 

°5 

DC  motor,  brush  commutation  loss 

kW 

E4 

E4 

DC  motor,  back-emf 

pu 

ES 

Eo 

DC  motor,  terminal  voltage 

pu 

ER 

E 

r 

AC  machine,  air  gap  voltage 

pu 

E5 

E5 

AC  machine,  excitation  voltage 

pu 

EX 

E 

X 

Controller,  commutation  voltage  drop 

V 

F. 

1 

Vehicle,  inertia  force 

lb 

F 

g 

Vehicle,  gravity  (grade)  force 

lb 

R FORCE i 

K 

Vehicle,  tire  friction  force 

lb 

F 

l w 

Vehicle,  aerodynamic  force  (inch  wind) 

lb 

FT 

F 

\fehicle,  required  tractive  force 

lb 

Fb 

Fo 

Tractive  force/ dc  motor  base  tract,  force 

pu 

GRADE 

G 

1 

I 

a, 

Grade  of  roadway  (profile) 

rad 

IAB 

DC  motor,  base  armature  current  (with  field  weak. ) 

U 

C 

DC  motor,  armature  current 

pu 

C-2 


(continued) 


Symbol  Symbol 

(Program)  (Equations) 

Definition 

Units 

IFR 

DC  motor,  shunt  field  current  with  field  weak 

pu 

IFDC 

r 

if 

DC  motor,  total  field  current 

pu 

IFACR 

ac 

AC  machine,  field  current  for  air  gap  voltage  (E^) 

pu 

14  B 

ac 

V 

AC  machine,  armature  current  after  rewind 

A 

14 

AC  machine,  armature  current  (per  phase) 

pu 

14  F 

s 

Fundamental  component  of  1^ 

pu 

15 

E 

AC  machine,  required  field  current 

pu 

KAGR 

K 

Slope  of  air  gap  line,  reduced  speed 

pu 

K2 

ag 

r 

KZ 

Flywheel,  moment  of  inertia  coefficient 

pu 

K3 

K3 

Flywheel,  maximum  speed  coefficient 

pu 

KFACT 

Kfact 

Ratio  of  slopes  of  air-gap  lines,  Kingsley  Method 

pu 

PINW 

P. 

in 

Charging  power,  ac  motor 

kW 

PRDC 

w 

p 

Loss  power,  dc  motor  rotor 

kW 

PSDC 

rdc 

P 

Loss  power,  dc  motor  stator 

kW 

PRAC 

Sdc 

P 

Loss  power,  ac  machine  rotor 

kW 

PSAC 

rac 

P 

Loss  power,  ac  machine  stator 

kW 

PDC  OUT 

s 

ac 

Pdc 

DC  motor,  rate  of  cooling  by  forced  air 

kW 

PAC  OUT 

out 

p 

AC  machine,  rate  of  cooling  by  forced  oil 

kW 

P 

ac  , 
out 

P 

Total  power  from  flywheel 

kW 

C-3 


(continued) 


Symbol  Symbol 


(Program  ) 

(Equat  ions) 

Definition 

Units 

P2 

P2 

Propulsion  power  at  wheels 

kW 

P3 

P3 

Total  losses  of  propulsion  system 

kW 

P4 

P4 

Transmission  losses 

kW 

P5 

P5 

Total  dc  motor  loss 

kW 

P6 

P6 

Total  controller  loss 

kW 

P7 

P7 

Total  ac  machine  loss 

kW 

PPDR 

P , 
pdr 

Power,  alternator  to  controller 

kW 

POW 

POW 

Total  alternator  power 

kW 

QPDR 

^pdr 

Fundamental  alternator  reactive  power 

kvar 

Q 

Q 

Total  alternator  fundamental  reactive  power 

kvar 

Q1 

Q, 

Reactive  power,  auxiliary  load 

kvar 

QH 

Qh 

Reactive  power,  alternator  - higher  freq.  comp. 

kvar 

RADC 

R 

a , 
dc 

DC  motor  , armature  circuit  res.  at  20°  C w / field  weak  ft 

R 1 

Ri 

Controller,  snubber  losses 

kW 

R2 

R2 

Controller,  SCR  losses 

kW 

S 

S 

Total  kVA  output  of  alternator 

kVA 

SI 

S1 

Fundamental  mode  of  total  alternator  output 

kVA 

SQ1 

SQ 

DC  motor,  armature  current  parameter 

— 

SLOPE B 

Slope  B 

AC  machine  air  gap  voltage/  corres.  field  current 

pu 
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(continued) 


Symbol 

Symbol 

(Program) 

(Equations) 

Definition 

Units 

TRDC 

T 

r 

ac 

Temperature  rise,  dc  motor  rotor 

°C 

TSDC 

T 

Sdc 

Temperature  rise,  dc  motor  stator 

°C 

TGDC 

T 

gdc 

Temperature  rise,  dc  motor  air  gap 

°c 

TRAC 

T 

r 

ac 

Temperature  rise,  ac  machine  rotor 

°c 

TSAC 

T 

s 

ac 

Temperature  rise,  ac  machine  stator 

°c 

TGAC 

T 

gac 

Temperature  rise,  ac  machine  cooling  oil 

°c 

T6 

T6 

Temperature  rise,  controller,  heat  sink  (fins) 

°c 

T6T 

T6 

t 

t 

Temperature  rise,  controller,  SCR  junction 

°c 

TIME 

Time  from  start  of  a run 

s 

DELTAT  At 

Computation  step  size 

s 

VMPH  - 

VMPHY 

A V 

Change  in  vehicle  velocity  over  At 

mi/  h 

V 

V 

Vehicle  velocity  (from  profile) 

ft/  s 

v<t> 

V 

o 

Vehicle  velocity/ vehicle  vel.  at  dc  motor  base 

speed 

pu 

H WIND 

V 

w 

Vdc 

Wind  speed  (headwind  is  positive) 
Terminal  voltage,  dc  motor 

mi/  h 
V 

VT 

vt 

Terminal  voltage,  ac  machine  (l -n) 

V 

VTB 

v,b 

Base  ac  machine  term,  voltage  (i  -n)  after  rewind 

V 

XDSAT 

xd  t 
sat 

Direct  axis  saturated  synchronous  reactance 

pu 
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Symbol 

Symbol 

(Program ) 

(Equations ) 

Definit  ion 

Units 

A 

a 

Firing  angle  for  controller 

deg 

D 

6 

Power  angle  for  ac  machine 

deg 

NGB 

gb 

Gearbox  efficiency 

deg 

NDC 

^dc 

DC  motor  efficiency 

NPD 

^pdr 

Controller  efficiency 

NAC 

^ac 

AC  machine  efficiency 

PHIDC 

dc 

Total  dc  flux  - conductor  product 

pu 

T8 

9 

AC  machine  power  factor  angle 

deg 

W<£ 

(J 

o 

Flywheel  and  alternator  speed/  maximum  speed 

pu 

** 
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charging  mode  data  : (Station  A,  Initial) 
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APPENDIX  E,  MANUFACTURERS'  DATA  FOR  MAJOR  COMPONENTS 


E - 1 


Part  I ALTERNATOR  DATA 


Source:  Westinghouse  Aerospace  Electrical  Division,  Lima,  Ohio 


15 


Generator  part  number:976J936- 1 0 


Service: 

Physical  Constants: 


Used  in  U.  S.  Navy  S3-A  aircraft 

Weight  42.  0 lb 

Size  8.  0 D x 9.  7 long 

Drawing  No.  946F907 

Rating  75  kVA  at  12,000  r/ min 

Number  of  poles  4 

Insulation  ML-rated  at  200° C 

for  l,000  h (continuous) 
Cooling  Spray  oil  at  4 gal/  min 

Characteristics  3 -phase,  brushless 

generator,  115  V(l-n) 


Weight  breakdown  in 

lb: 

Stator  steel 

11.0 

Rotor  steel 

6.  9 

Rotor  copper 

2.  9 

Stator  copper 

4.  1 

Frame  plus  insulation 

17.  1 

Total  weight 


42.  0 


E-2 


(continued) 


Symbol 

PU 

Ohms 

Armature  resistance  (200° C) 

! 

a 

0.  0495 

0.  0268 

Direct  synchronous  reactance  (unsaturatedy  X , 

d 

3.  4 

1.  84 

Quadrature  synchronous  reactance 

X 

q 

1.  24 

0.  67 

Direct  subtransient  reactance 

i t 

xd 

0.  145 

0.  078 

Armature  leakage  flux  reactance 

X 

a 

0.  085 

0.  046 

Losses  in  Watts 


No 

Load 

1 00%  Load 
(75.  8 kW) 

Loss  Model 
Term 

Stator  copper  loss 

0 

3938 

A. 

Rotating  field  copper 

251 

1687 

A2 

Exciter  loss 

260 

300 

A2 

Stator  core  loss 

1230 

1255 

A3 

Pole  face  loss 

325 

346 

A3 

Windage  loss 

2045 

2045 

\ 

Stray  load  loss 

207 

753 

A5 

Fan  loss 

— 

— 

— 

Total  431  8 10324  P? 

75, 800 


Efficiency  = 


75, 800  + 10, 324 


x 1 00 


88.  0%  at  rated  load 


E-3 


Part  II  DC  TRACTION  MOTOR  DATA 


Source:  ASEA  manual  for  70  hp  dc  motor  LAU203M 


Rating 

7o  HP 

Base  speed 

273>o  Rpm 

Armature  voltage 

42o  Volts 

Armature  current 

14o  Amps 

Field  voltage 

24  volts 

Field  current 

12  Amps 

Number  of  poles 

4 

Enclosure 

Drip-proof 

Standards 

Nema 

Insulation  class 

F 

The  motor  is  an  ordinary,  ventilated  version  with  two  free  shaft  ends 
It  has  been  provided  with  specially  strengthened  feet  suitable  for 
ceiling  mounting. 

The  internal  fan  is  mounted  opposite  the  commutator,  and  the  cool- 
ing air  is  drawn  in  radially  at  the  commutator  end  and  is  exhaust- 
ed radially  at  the  fan  end.  Furthermore  it  is  provided  with  a 
special  air  duct  at  the  commutator  end  for  tube  connection  to  a 
separate  blower. 

The  motor  is  fed  from  a 3-phase,  39oV,  6oHz  supply  through  a fully 
controlled  thyristor  controller.  The  shunt  excitation  is  taken  from 
a 24  Volt  DC-supply.  This  supply  is  controlled  in  order  to  increase 
the  motor  speed  to  317o  rpm. 


t 


E-  4 


Cross  Section  of  70HP-DC  Traction  Motor 
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Table  E-l 

Parts  List  For  70  HP- 
DC  Traction  Motor.  (See  p.  E-4) 
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Table  E - 2_ 

Test  Data  For  70  HP-DC  Traction  Motor 


Cw.io-.r  Ait  A Inc. 

A New.  King  St. 
White  Plains 
N.  Y.  10602 
U.  S.  A. 

Up* 

LAU  203M 

No 

2184423 

Output 

70  h.p. . kw 

Service 

cont . 

Voltoge 

420  v 

Current 

140  A 

Kind  of  mocMn* 

ComD.  motor 

Excitation 

24  v 

Speed 

2730  rpm 

Vo  1 'a  j • Cvrfnt  A Input 

V orm.  | e«c.  kW 

Speed  Co»pof  Efficiency 

rpm  k'«V  X 

No-lood 


420 

420 


5,1 

5,1 


12,5 

12,5 


3000 

■3000 


clockwise 
counterclcck . 


load 


420  140 

420  140 


12,5 

12,5 


2680 

2680 


clockwise 

counterclock. 


Re«i«tan;e 

between 

terminal, 

ohm. 

Winding 

a»  20 

•c 

Arm. 

0,0646 

0,0327 

Comm. 

0,0111 

1,67 

Serie, 

Shunt 

Insul.  resistance  to 
frame,  megohms 

100 

Temperature 

rise 

•c 

Winding, 
port  of  machine 

After  continuous  test  at  26S0  rpm 

420  V 140  A 12,5  A exc. 

After  at  rmp 

V A A exc. 

Meotured  by 
thermometer 

B 

In  reiUtance 

embedded 

Meosured  by 
thermometer 

By  Inereote 
In  reilitonce 

8 r embedded 
thermometer 

Arm. 

69 

37 

Comm. 

29 

38 

Serl«, 

Shunt 

Commutator 

32 

High  voltage  te»t. 

LtC.  winding  to  frame 

A/m.  and  cone,  winding  to  frome 

lefween  exc.  winding 

2000  v 1 

min. 

2000  V 1 min. 

2000  V 

l min- 

Over-  role 
tests 

E»c*u  voltage 

V ( 'L  incr.) 

min. 

uuu  current 

A ( y.  incr.)  min. 

£>»•'»  peed 

rpm  ( 

U incr.)  min. 

I 
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o 


u 

rt 


CNJ 

I 

H 

W 

•H 

PU 


ex 


<D 

> 

X 

3 

U 


0) 

u 

C 

3 

6 

X 

o 

<4-1 

X 

0) 

Cp 
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WINDING  DATA 
ARMATURE 

43  slots,  129  commutator  bars 

43  coils  (l-f-2-f-l)  turns  of  wire  3 par.  d=2,  o(insulated  with  a double 
coating  of  enamel). 

Weight  of  copper  : 24  lbs. 

SHUNT  COIL 

4 coils,  15o  turns  of  v/ire  d=  2,24  (insulated  with  a double  coating 
of  enamel). 

Weight  of  copper:  3o  lbs.  ■ ' 

Resistance  : 1,7  ohms  cold. 

SERIES  COIL 

4 coils,  12  turns  of  wire  3 x 5 m (Varnish  bonded  glass  covered 
rectangular  copper). 

Weight  of  copper:  11  lbs. 

INTERPOLE  COIL 

4 coils,  51  turns  of  v/ire  3 x 5 mm  (Varnish  bended  glass  covered 
rectangular  copper). 

Weight  of  copper  : 33  lbs. 
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Table  E-3 

Power  Pad  Thyristor  Data 
SOLID-STATE  CONTROLLER  DATA 


Source : 


Power  Semiconductors, 


Inc.  publication: 


9 PSI-33(909) 


TYPICAL  CONTINUOUS  USABLE  D.  C.  RATINGS 


1500  LFM 

1 DC 

KW  Output  (1  DC  Leg  X3  X V OC) 

DC  MOTOR  DRIVE  - 1.15  SERVICE  FACTOR 

Fan  Cooled 

Per  Leg 

@ 300  V DC 

@600  V DC 

@900  V DC 

@ 240  V DC 

@500  V DC 

@ 700  V DC 

@900  V DC 

F-180 

125 

110 

225 

335 

85 

200 

285 

365 

F-220 

140 

125 

250 

375 

110 

225 

315 

405 

F-300 

170 

150 

305 

460 

130 

275 

385 

500 

F-400 

210 

190 

375 

565 

160 

340 

475 

610 

F-500 

230 

205 

415 

620 

180 

370 

520 

670 

F-600 

250 

225 

— 

— 

195 

— 

— 

— 

Based  on:  6"  PSI -3007  Aluminum  heatsink  (0.08°  c/w  H.  S.  - Amb.  Thermal  Impedance) 

RATINGS: 

TYPE  PSI 

E or  F-180 

E or  F-220 

E or  F 300 

Eot  F-400 

E or  F-500 

E or  F-600 

RMS  Fwd.  Cur. 

Irms 

180A 

220A 

300A 

400A 

500A 

600A 

AVE  Fwd.  Cur.  (180°  Cond.)  Iave 

115A 

140A 

190A 

250A 

320A 

375A 

Fwd.  Volt  Drop  @ 500A  Peak  Vfm 

3.00V 

2.60V 

2.00V 

1.50V 

1.35V 

1.20V 

l2t  for  fusing  @ 8.3  mS 

i2t 

35.000 

35.000 

65,000 

80.000 

100.000 

125,000 

Surge  Cur.  (1)  cycle 

Ifm 

3000A 

3000A 

4000A 

4500A 

5000A 

5500A 

Table  E-3 


E - 1 1/E  - 1 2 


(Continued) 


100mA(typ.)  'M 


*710 


•For  Special  Selections  consult  Factory  - 
dv/dt  up  to  5,000  V/u  S available. 


TYPE  E or  F-600 

E-180,  220,  300, 400  or  500 
F-180,  220,300,  400  or  500 

-5 

-5 

-5 

-6 

-6 

-6 

•0  OO  00 
1 1 1 

-10 

-10 

HIGHER 

-12 

-12 

'JOT  AVA 
-14 
-14 

ILABLE 

N. 

-16 

A. 

-18 

PEAK  FWD.  & REV.  VOLTS 
-40  + 0 + 125°  c 

Vfom 

Vrom 

500 

600 

800 

1000 

1200 

1400 

1600 

1800 

PEAK  NON.  REP.  TRANSIENT  VOLT. 

600 

700 

900 

1100 

1300 

1500 

1700 

1900 

PEAK  FWD.  & REV.  LEAKAGE 
CURRENT- MILLIAMPERE 

Ifom 

Irom 

15 

15 

15 

10 

10 

10 

10 

10 

GATE 
Peak  Fwd.  Current 


Peak  Power,  Watts 


Ave  Power,  Watts 


Current  to  Fire,  25°  C 


Volts  to  Fire,  25°  C 


Volt,  not  to  Fire,  125°  C 


THERMAL  IMP.  J-C 


THERMAL  IMP.  C-HS 


TURN  ON  TIME,  If  = 5A 


TURN  OFF  TIME  @ 125°  C 


MIN.  dv/dt  @125°  CLIN. 


MIN.  di/dt  @ 125°  C 


HOLDING  CURRENT,  25°  C 


LATCHING  CURRENT,  25°  C 


MOUNTING  FORCE  ±10% 


Igfm 


Pgm 


12W 


Pg(av) 


2.5W 


150mA 


TRANSIENT  THERMAL  IMPEDANCE 


J unction  to  Case  Double  S_i&edL  .Coole dj| 


Vgt 


Vgd 


0.15V 


0.09°  c/w 


0c-s 


0.03O  c/w 


T on 


10  uS 
(typical) 


Toff 


150  uS 
(typical) 


dv/dt 


100V/oS* 


di/dt 


50  A/uS* 


500mA 


1200  lbs. 


mm 


Kmaximum  forward  voltage  drop  @ 2 5 cm 


APPENDIX  F,  REPORT  OF  INVENTIONS 


The  objective  of  this  work  was  to  develop  a computer 
simulation  of  electric  propulsion  systems  for  flywheel  energy 
storage  vehicles.  The  components  that  were  used  to  characterize 
the  propulsion  system  were  standard  components,  therefore 
the  work  resulted  in  no  discoveries  or  inventions. 
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